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Introduction 



Since the late 1940s, nanometer-sized magnetic particles have generated continuous 
interest as the study of their properties has proved to be scientifically and technolog- 
ically very challenging. In particular it was recognized that the ferromagnetic state, 
with a given orientation of the particle moment, has a remanent magnetization if the 
particle is small enough. This was the starting point of huge permanent magnets and 
magnetic recording industries. However, despite intense activity during the last few 
decades, the difficulties in making nanoparticles of good enough quality has slowed 
the advancement of this field. As a consequence, for 50 years, these applications 
concentrated above and then near the micrometer scale. In the last decade, this has 
no longer been the case because of the emergence of new fabrication techniques 
which have led to the possibility of making small objects with the required struc- 
tural and chemical qualities. In order to study these objects new techniques were 
developed such as magnetic force microscopy, magnetometry based on micro-Hall 
probes or micro-SQUIDs. This led to a new understanding of the magnetic behavior 
of nanoparticles, which is now very important for the development of new fundamen- 
tal theories of magnetism and in modehng new magnetic materials for permanent 
magnets or high density recording. 

In order to put this review into perspective, let us consider Fig. I.l which presents 
a scale of size ranging from macroscopic down to nanoscopic sizes. The unit of 
this scale is the number of magnetic moments in a magnetic system. At macroscopic 
sizes, a magnetic system is described by magnetic domains (Weiss 1 907) [ 1] which are 
separated by domain walls, magnetization reversal occurs via nucleation, propagation 
and annihilation of domain walls (see the hysteresis loop on the left in Fig. I.l which 
was measured on an individual elUptic CoZr particle of 1 /xm x 0.8 fim, and a 
thickness of 50 nm). Shape and width of domain walls depend on the material of 
the magnetic system, on its size, shape and surface, and on its temperature [2]. The 
material dependence of the domain walls has motivated the definition of two length 
scales: (i) the domain wall width S defined hy S = ^ A/K and (ii) the exchange 
length A defined by A = vG4/Ms where A is the exchange energy, K is the crystalline 
anisotropy constant and Ms is the spontaneous magnetization. Qualitatively, the first 
definition shows that anisotropy energy favors a thin wall, while the exchange energy 
favors a thick wall. For very small crystalline anisotropy, the first definition suggests 
an infinite domain wall width which has a large total energy. This is due to the 
magnetostatic energy term which can be reduced by subdividing the ferromagnetic 
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Fig. I.l Scale of size which goes from macroscopic down to nanoscopic sizes. The unit of 
this scale is the number of magnetic moments in a magnetic system (roughly corresponding to 
the number of atoms). The hysteresis loops are typical examples of magnetization reversal via 
nucleation, propagation and annihilation of domain walls {left), via uniform rotation {middle), 
and quantum tunneling {right). 
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crystal into domains. Therefore, for very small crystalline anisotropy, the domain 
wall width is of the order of magnitude of the exchange length A. Both length scales 
can range from sub-micrometer scales in alloys to atomic scales in rare earth systems. 

When the system size is of the order of magnitude of 6 or A, the formation of domain 
walls requires too much energy. Therefore, the magnetization remains in the so- 
called single domain state. ^ Hence, the magnetization might reverse by uniform 
rotation, curling or other non uniform modes (see hysteresis loop in the middle of 
Fig. I.l). In this review, we discuss mainly this size range where the physics is rather 
simple (Sects. 2 and 3). 

For system sizes well below S and A, one must take into account explicitly the mag- 
netic moments (spins) and their couplings. The theoretical description is complicated 
by the particle's boundaries [3, 4, 5]. 

At the smallest size (below which one must consider individual atoms and spins) 
there are either free clusters made of several atoms [6, 7] or molecular clusters which 
are macro- molecules with a central complex containing magnetic atoms. In the last 
case, measurements on the Mni2 acetate and Fcg molecular clusters showed that the 
physics can be described by a collective moment of spin S=\0 (Sect. 4.1). By means 
of simple hysteresis loop measurements, the quantum character of these molecules 
showed up in well defined steps which are due to resonance quantum tunneUng 
between energy levels (see hysteresis loop on the right in Fig. I.l). 
In the following sections, we review the most important theories and experimental 
results concerning the magnetization reversal of single-domain particles and clusters. 
Special emphasis is laid on single particle measurements avoiding compUcations due 
to distributions of particle size, shape etc. Measurements on particle assemblies has 
been reviewed in |8]. We mainly discuss the low temperature regime in order to 
avoid spin excitations. 

In Sect. 1 , we briefly review the commonly used measuring techniques. Among them, 
electrical transport measurements. Hall probes and micro-SQUID techniques seem 
to be the most convenient techniques for low temperature measurements. Sect. 2 
discusses the mechanisms of magnetization reversal in single domain particles at 
zero kelvin. The influence of temperature on the magnetization reversal is reported 
in Sect. 3. Finally, Sect. 4 shows that for very small systems or very low temperature, 
magnetization can reverse via tunneling. 



In the theory of micromagnetism, the single domain state describes the state where the magnetization is 
perfectly ahgned [2] whereas experimentalists mean often a state without domain wall. 



1 Single particle measurement 
techniques 

The following sections review commonly used single particle measuring techniques 
avoiding complications due to distributions of particle size, shape etc. which are 
always present in particle assemblies [8]. Special emphasis is laid on the micro- 
SQUID technique and the developed methods which allowed the most detailed studies 
at low temperatures. 



1.1 OVERVIEW OF SINGLE PARTICLE MEASUREMENT 
TECHNIQUES 

The dream of measuring the magnetization reversal of an individual magnetic par- 
ticle goes back to the pioneering work of Neel [9, 10]. The first reahzation was 
published by Morrish and Yu in 1956 [11]. These authors employed a quartz- fiber 
torsion balance to perform magnetic measurements on individual micrometer sized 
7— Fe203 particles. With their technique, they wanted to avoid the complication of 
particle assembhes which are due to different orientations of the particle's easy axis 
of magnetization and particle-particle dipolar interaction. They aimed to show the 
existence of a single-domain state in a magnetic particle. Later on, other groups 
tried to study single particles but the experimental precision did not allow a detailed 
study. A first breakthrough came via the work of Knowles [12] who developed a 
simple optical method for measuring the switching field, defined as the minimum 
applied field required to reverse the magnetization of a particle. However, the work of 
Knowles failed to provide quantitative information on well defined particles. More 
recently, insights into the magnetic properties of individual and isolated particles 
were obtained with the help of electron holography [13], vibrating reed magnetome- 
try [14], Lorentz microscopy [15, 16], magneto-optical Kerr effect [17], and magnetic 
force microscopy [18, 19]. Recently, magnetic nano-structures have been studied 
by the technique of magnetic hnear dichroism in the angular distribution of photo- 
electrons or by photoemission electron microscopy [20, 21]. In addition to magnetic 
domain observations, element- specific information is available via the characteristic 
absorption levels or threshold photoemission Among all mentioned techniques, 
most of the studies have been carried out using magnetic force microscopy at room 



'We refer to the literature concerning other domain observation techniques [1]. 
4 
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temperature. This technique has an excellent spatial resolution but time dependent 
measurements are difficult due to the sample-tip interaction. 
Only a few groups were able to study the magnetization reversal of individual nanopar- 
ticles or nanowires at low temperatures. The first magnetization measurements of 
individual single-domain nanoparticles and nanowires at very low temperatures were 
presented by Wernsdorfer et al. [22]. The detector (a Nb micro-bridge- DC-SQUID) 
and the studied particles were fabricated using electron-beam hthography. Coppinger 
et al. [23J investigated the magnetic properties of nanoparticles by resistance mea- 
surements. They observed the two-level fluctuations in the conductance of a sample 
containing self-organizing ErAs quantum wires and dots in a semi-insulating GaAs 
matrix. By measuring the electrical resistance of isolated Ni wires with diameters 
between 20 and 40 nm, Giordano and Hong studied the motion of magnetic domain 
waUs [24, 25]. Other low temperature techniques which may be adapted to single 
particle measurements are Hall probe magnetometry [26, 27, 28J, magnetometry 
based on magnetoresistance [29, 30, 31] or spin-dependent tuimeling with Coulomb 
blockade [32, 33]. At the time of writing, the micro-SQUlD technique allows 
the most detailed study of the magnetization reversal of nanometer-sized particles 
[34, 35, 36, 37, 38, 39]. The following section reviews the basic ideas of the micro- 
SQUID technique. 



1.2 MICRO-SQUID MAGNETOMETRY 



The Superconducting Quantum Interference Device (SQUID) has been used very 
successfully for magnetometry and voltage or current measurements in the fields of 
medicine, metrology and science [40, 41]. SQUIDs are mostly fabricated from a 
Nb— AlOx— Nb trilayer, several hundreds of nanometers thick. The two Josephson 
junctions are planar tunnel junctions with an area of at least 0.5 /im^. In order to 
avoid flux pinning in the superconducting film the SQUID is placed in a magnetically 
shielded environment. The sample's flux is transferred via a superconducting pick 
up coil to the input coil of the SQUID. Such a device is widely used as the signal 
can be measured by simple lock-in techniques. However, this kind of SQUID is not 
well suited for measuring the magnetization of single submicron- sized samples as 
the separation of SQUID and pickup coil leads to a relatively small coupling factor. 
A much better couphng factor can be achieved by couphng the sample directly with 
the SQUID loop. In this arrangement, the main difficulty arises from the fact that 
the magnetic field applied to the sample is also applied to the SQUID. The lack of 
sensitivity to a high field apphed in the SQUID plane led us to the development of 
the micro-bridge-DC-SQUID technique [22] which allows us to apply several Tesla 
in the plane of the SQUID without dramatically reducing the SQUID's sensitivity. 
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1.2.1 Choice of SQUID configuration 

The main criteria for the choice of the micro-SQUID configuration were an easy 
couphng to a mesoscopic sample, a simple fabrication, a simple mode of operation, 
robustness and stability, the desired temperature range, and operation in high magnetic 
fields (in particular for fields applied in the SQUID plane). These criteria led to the 
use of micro bridge junctions instead of the commonly used tunnel junctions. 
The Josephson effect in micro bridge junctions has first been suggested in 1964 by 
Anderson and Dayem [42]. These superconducting weak links seemed to be very 
promising in order to design planar DC-SQUIDs with a one step thin film technology. 
However, Dayem bridges exhibit a Josephson current-phase relation only when their 
dimensions are small compared to the coherence length ^. Such dimension were 
difficult to reach at those days. Nowadays, electron beam hthography allows one to 
directly fabricate reliable micro-bridge Josephson Junctions made of materials Uke 
Al, Nb and Pb. 

1.2.2 Fabrication of micro-SQUIDs 

By using electron beam lithography, planar Al or Nb micro-bridge-DC-SQUIDs (of 
0.5 to 4 /xm in diameter) can be constructed (Fig. 1.1) [43, 44, 45]. Al SQUIDs can 
be obtained by evaporating a 20 to 50 nm thin film onto a PMMA mask, followed 
by standard lift-off techniques. The fabrication of Nb SQUIDs revealed to be more 
difficult. The direct evaporation of Nb onto a PMMA mask led to Nb SQUIDs 
of poor quality which manifested itself in a very low critical temperature of the 
superconductivity. Better Nb qualities were achieved by using UHV facilities and 
substrate temperatures of about 800 K, or higher, in order to grow 20 nm thin Nb 
films on Si substrates. ^ The Nb films were covered by 5 nm of Si. The Nb SQUIDs 
were then patterned with Reactive Ion Etching (RIE) using an AlO mask made by a 
electron beam lithography via a standard PMMA technique. 
An alternative method to fabricate micro-bridge-DC-SQUIDs has recently been pro- 
posed by Bouchiat et al. [46]. The new method is based on local anodization of 
3 to 6.5 nm thick Nb strip lines under the voltage biased tip of an Atomic Force 
Microscope. Micro-bridges junctions and SQUID loops were obtained either by 
partial or total oxidation of the Nb layer The first fabricated devices had about the 
same performance than micro-SQUIDs fabricated by electron beam lithography. The 
AFM-made SQUIDs should offer new features such as the fabrication at a chosen 
position allowing an optimized coupling to magnetic signals. In addition, we expect 
an increased intrinsic sensitivity: In the case of small magnetic clusters which are 
placed very close to the micro-bridge junctions, an improvement of the sensitivity 



^In order to have a SQUID which can be exposed to a high field appUed in the SQUID plane, the SQUID 
was made out of a very thin layer preventing flux trapping. In most cases, we used mostly 20 nm thick 
Nb layers allowing measurements of hysteresis loops in magnetic fields of up to 2 T. Such SQUIDs might 
work at fields higher than 10 T when using very thin (< 10 nm) Nb layers of high quality. 
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Fig. 1.1 Scanning electron micrograph of a Nb micro-bridge-DC-SQUID fabricated by elec- 
tron beam lithography. A Ni wire of diameter of about 90 nm was deposit on the SQUID 
(Sect. 2.2.2). 



of one to two orders of magnitude might be achieved due to the reduction of the 
micro-bridges size. It might allow us to detect the spin flips of about 100 magnetic 
moments. 

1.2.3 Magnetization measurements via critical current measurements 

The micro-bridge-DC-SQUID have a hysteretic V — I curve (Fig. 1.2): Ramping 
the current up from zero, the SQUID transits from the superconducting to the normal 
state at the critical current 4. Due to Joule heating in the normal state, the SQUID 
stays resistive down to currents much smaller than Jc. This hysteretic V — I curve 
made it impossible to use standard SQUID electronics or lock- in amplifier to read out 
the SQUID. Therefore, Benoit et al. developed a method consisting in measuring the 
critical current of the micro-SQUID [43, 44]. A computer-controlled circuit triggers 
simultaneously a current ramp (Fig. 1.3) and a 40-MHz quartz clock. As soon as a 
dV/ dt pulse is detected at the SQUID due to the transition from the superconducting 
to the normal state ^, the clock stops and the current is set to zero. The clock reading 
is transferred to the computer, and the cycle begins again. The critical current is 
proportional to the duration of the current ramp. The repetition rate is about 10 kHz, 
Umited by the time needed to settle the current. As the critical current Ic is a periodic 



^It is important to mention that the dV/dt pulse can be detected directly on the current biasing lead of the 
SQUID, i.e. for most cases it is sufficient to connect the SQUID with a single wire and the mass of the 
cryostat. 
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Fig. 1.2 Oscilloscope reading of a voltage versus current curve of a Nb micro-bridge-DC- 
SQUID. The SQUID transits from the superconducting to the normal state at the critical current 




Fig. 1.3 Current injected into SQUID loop. First, the current is increased quickly up to a 
current /c — e which is close to the critical current Ic- Then, the current is ramped up to 7c- As 
soon as a dV/ dt pulse is detected at the SQUID due to the transition from the superconducting 
to the normal state the current is set to zero. 1^ is proportional to the duration of the current 
ramp At. The repetition rate of the cycle is up to 10 kHz. 
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Fig. 1.4 Critical current Ic of a Al micro-SQUID as a function of a magnetic field applied 
perpendicular to ttie SQUID plane. For the feedback mode an external field perpendicular to 
the SQUID loop is applied which keeps the critical current constant at a working point, for 
example point P or P'. For the cold mode method the SQUID is biased close to the critical 
current so that it is in state A or B, respectively of positive or negative flux jumps induced by 
the magnetization reversal. 



function of the flux going through the SQUID loop (Fig. 1.4), one can easily deduce 
the flux change in the SQUID loop by measuring the critical current. The sensitivity 
achieved by the critical current measurement technique was about 10~^^o / VHz for 
Al SQUIDs and IQ-^^o/v^ for Nb SQUIDs ($o = h/2e = 2 x 10"^^ Wb). 
In order to have good magnetic flux coupling, the mesoscopic systems, for in- 
stant the magnetic particles, is directly placed on the SQUID loop (Fig. 1.5) [45]. 
The SQUID detects the flux through its loop produced by the sample's magnetiza- 
tion. For hysteresis loop measurements, the external field is applied in the plane 
of the SQUID (Fig. 1.5), thus the SQUID is only sensitive to the flux induced 
by the stray field of the sample's magnetization. The flux sensitivity of the criti- 
cal current mode allowed us to detect the magnetization reversals corresponding to 
Hz{l0^^^^emu/\/Hz), i.e. the magnetic moment of a Co nanoparticle 
with a diameter of 5 nm. The time resolution was given by the time between two 
measurements of the critical current. In this case, the achieved time resolution was 
100 /is. 

1.2.4 Feedback mode for hysteresis loop measurements 

For the detection of flux variations larger than $o /2, the direct critical current method 
(Sect. 1.2.3) is complicated by the fact that Ic is a non monotonous function of the 
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Fig. 1.5 Schematic drawing of a planar micro-bridge-DC-SQUID on whicti a ferromagnetic 
particle is placed. The SQUID detects the flux through its loop produced by the sample 
magnetization. Due to the close proximity between sample and SQUID a very efficient and 
direct flux coupling is achieved. 



flux (Fig. 1 .4) which goes through the SQUID loop. In this cases, we used a feedback 
mode consisting in applying an external field perpendicular to the SQUID loop 
which compensates the stray field variations keeping the critical current constant 
at a working point as shown in Fig. 1.4. After each Ic measurement, the computer 
calculates the difference between a working point and Ic, multiplies it with a feedback 
factor, and adds it to the feedback (or compensation) field. The feedback field gives 
then the stray field variation of the sample's magnetization. It can easily be given 
in unites of $o- Knowing the couphng factor between particle and SQUID loop, 
absolute magnetic moment measurements are possible. 

The time resolution of the feedback mode depends on the rate of the 1^ measurements. 
For rates of about 10 kHz, we are able to follow flux variations in the milUsecond 
range. 



1.2.5 Cold mode method for magnetization switching measurements 

For studying the magnetization switching of nanoparticles, we developed a special 
mode, called the cold mode method [36, 37, 38, 39], which are much faster and more 
sensitive than the previous presented modes. The achieved flux sensitivity allowed 
us to detect the magnetization switching corresponding to 10^/iB(10~^'^emu), i.e. 
the switching of a Co nanoparticle with a diameter of 2 to 3 nm. The time resolution 
of the switching detection reached the nanosecond range. 

The cold mode method is also specially adapted for studying the temperature depen- 
dence of the magnetization reversal and macroscopic quantum tunneling of magne- 
tization. Indeed, the main difficulty associated with the SQUID detection technique 
lies in the Joule heating when the critical current is reached. After the normal state 
transition at the critical current, the SQUID dissipates for about 100 ns which slightly 
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heats the magnetic particle coupled to the SQUID. This problem can be solved by 
the cold mode method which uses the SQUID only as a trigger [36, 37, 38, 39]. 
The cold mode consists in biasing the SQUID close to the critical current while a field 
is applied perpendicular to the SQUID plane so that the SQUID is in state A (B) for 
a positive (negative) flux jump (Fig. 1 .4) which can be induced by the magnetization 
reversal of a particle coupled to the SQUID loop. The magnetization reversal then 
triggers a transition of the SQUID from the superconducting to the normal state, 
i.e. a dV/dt pulse can be detected on the current lead biasing the SQUID. Our 
SQUID electronics allows us to detect the dV/ dt pulse some nanoseconds after the 
magnetization reversal, allowing very precise switching field measurements. Because 
this method only heats the sample after the magnetization reversal, we called it the 
cold mode method. 

Another advantage of this mode is that the sample does not interfere with the rf-noise 
which is induced in oxide layer Josephson junctions of conventional SQUIDs because 
the hysteretic micro-SQUID is in the superconducting state before the magnetization 
reversal. 

Finally, the cold mode method is very important for studying macroscopic quantum 
tunneUng of magnetization. Quantum theory predicts that the escape rate from a 
metastable potential well by quantum tunneling is strongly reduced by the coupling 
of the magnetic system with its environment. Therefore, the measuring device must 
be weakly coupled to the magnetic particle. However, in order to measure the 
magnetization reversal, the SQUID must be strongly coupled to the magnetic particle 
which hinders the possibility of quantum tunnehng. This problem can be solved by 
using the cold mode method. In order to show this schematically. Fig. 1.6 represents 
two energy potentials, one is the double well potential of the particle, the other is 
the periodic potential of the SQUID. Before the magnetization reversal, both systems 
are in a metastable state: the particle because of an applied field which is close to 
the switching field and the SQUID because of a current through the SQUID loop 
which is close to the critical current. When the particle overcomes the saddle point or 
tunneUng through the energy barrier, its magnetization rotates by only few degrees. 
For this starting process of the magnetization reversal, the coupling between particle 
and SQUID can be arranged to be very small. Afterwards, the particle falls into 
the lower well which impUes a rotation of magnetization of up to 180°. During 
this process, the coupling between particle and SQUID should be strong enough to 

'^In order to illustrate the dipolar couplings, let us consider the energy scales involved. For most of the 
particles measured so far below 1 K, the energy barrier height from the metastable state up to the saddle 
point is of the order of a few Kelvins whereas that from the lower state up to the saddle point is between 
10^ and 10® K. These energy scales should be compared with the energy necessary to drive the SQUID 
out of its metastable superconducting state which is of the order of a few Kelvin. Therefore, only a small 
energy transfer is necessary to measure the magnetization reversal. In addition, a proper orientation of 
the easy axis of magnetization with respect to the SQUID loop can further reduce the coupling during the 
first stage of the magnetization reversal. In the case of an easy axis of magnetization perpendicular to the 
current direction in the SQUID wire (Fig. 1.5), the coupling factor between SQUID loop and particle is 
proportional to (1 - cos f ), where f is the angle between the direction of magnetization and its easy axis. 
Therefore, the coupling is very weak at the first stage of magnetization reversal where ip is small. 
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Fig. 1.6 Energy scheme of the cold mode method. After the particle overcomes the saddle 
point or tunnels through the energy barrier, it falls into the lower potential well releasing 
energy. A very small fraction AE of this energy is transferred to the SQUID and drives the 
SQUID out of its metastable superconducting state. 



drive the SQUID out of its metastable state. The corresponding transition from the 
superconducting into the normal state is easily measurable for a hysteretic SQUID. 
The main disadvantage of the cold mode is that only the switching field of magne- 
tization reversal can be measured and not the magnetization before and during the 
magnetization reversal. 

1.2.6 Blind mode method for three dimensional switching field measurements 

A disadvantage of the micro-SQUID technique is that it does not function properly 
when a significant field is applied perpendicular to the SQUID loop. It works 
also only below the superconducting critical temperature of Nb (Tc « 7 K for our 
SQUIDs). These facts limited us to 2D measurements and to T < Tc in the first 
experiments [35, 36, 37]. However, we showed recently that full-three-dimensional 
measurements can be done by using an indirect method [38]. In addition, this 
technique allows us to study the magnetization reversal for T > Tc- 
Let us consider Fig. 1.7 showing switching field measurements (Sect. 2) for in- 
the-SQUID-plane applied fields, detected using the cold mode (Sect. 1.2.5). The 
three dimensional switching field measurements and the studies as a function of 
temperature can be done using a three step method which we call the blind mode 
(Fig. 1.7): 



1 . Saturation: The magnetization of the particles is saturated in a given direction 
(atT = 35mK). 
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2. Testing: A test-field is applied at a temperature between 35 mK and 30 K ^ 
which may or may not cause a magnetization switching. 

3. Probing: After cooUng to T = 35 mK, the SQUID is switched on and a field 
is swept in the plane of the SQUID to probe the resulting magnetization state 
using direct critical current measurements, the feedback mode, or the cold 
mode. 



If the SQUID detects a magnetization jump in step (3), this means that the previously 
apphed test-field was weaker than the switching field for the direction being probed 
in step (2). The next interation will then be done with a stronger test-field. On 
the other hand, if the SQUID does not detect any magnetization jump in step (3), 
this means that the reversal occurred during step (2). The next interation will then 
be done with a weaker field. When choosing the new test-field with the help of a 
bisection algorithm, we needed about 8 repetitions of the three steps in order to get 
the switching field with good precision. This method allows us to scan the entire 
field space. * 

1.2.7 Micro-SQUID arrays 

We also use arrays of micro-SQUIDs as a magnetometer for macroscopic samples. 
There are three applications which were particular interesting: crystals of magnetic 
molecular clusters [47] (Sect. 4.1), nucleation and depinning of magnetic domain 
walls in thin films [48], and arrays of magnetic dots [49]. The procedure consists 
in placing a sample on top of an array of micro-SQUIDs so that some SQUIDs are 
directly under the sample, some SQUIDs are at the border of the sample, and some 
SQUIDs are beside the sample (Fig. 1.8). When a SQUID is very close to the sample, 
it is sensing locally the magnetization reversal whereas when the SQUID is far away, 
it integrates over a bigger sample volume. Therefore, depending on the sample, 
one can obtain more insight in the magnetization reversal than with conventional 
techniques which measure only the total magnetization. 

Our magnetometer works in the temperature range between 0.035 and 6 K and in 
fields up to 5 T with sweeping rates as high as 30 T/s, and a field stabiUty better 
than a microtesla. The time resolution is about 1 ms (Sect. 1.2.4) allowing short-time 
measurements. The field can be applied in any direction of the micro-SQUID plane 
with a precision much better than 0.1° by separately driving three orthogonal coils 
[45]. 



'Our highest temperatures of 30 K was only hrtiited by the cooling time. Below 30 K, we achieved cooling 
rates of few Kelvins per second. 

*It is worth mentioning that special precautions are necessary for anisotropics which are more complex 
than uniaxial (Sect. 2). 
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Fig. 1.7 Schematics for the blind mode method presenting the angular dependence of the 
switching field near the easy axis of a 3 nm Co cluster (full circles) (Sect. 2.1.2). First, the 
magnetization of the particle is saturated in a given direction (at T = 35 mK). Then, a test- 
field is applied, at a temperature between 35 mK and 30 K, which may or may not cause a 
magnetization switching. Finally, the SQUID is switched on (at T = 35 mK) and a field is 
applied in the plane of the SQUID to probe the resulting magnetization state. 



sam pie 




Fig. 1.8 Schematic representation of our magnetometer which is an array of micro- SQUIDs. 
Its high sensitivity allows us to study single crystals of the order of 10 to 500 /im which are 
placed directly on the array. 
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1.2.8 Scanning SQUID microscope 

The micro-SQUID technique has recently been used to build a scanning SQUID 
microscope [50]. The SQUID is designed by electron beam lithography at the apex 
of a silicon cantilever. The lever is attached to a force sensor, allowing to image 
magnetically, as well as topographically with a spatial resolution of 50 nm and a flux 
resolution of about 10~'^(j)o. The first application of this technique concerned the 
imaging of vortices in artificial networks. 

1.2.9 Outlook 

After the development of micro-SQUID technique in the early 90's |43, 44], the 
study of magnetization reversal in magnetic nano-structures began in 1993 [45]. 
The first studied systems were micrometer sized particles containing about 10^° 
magnetic moments. During the following years, the micro-SQUID technique has 
been improved to study smaller and smaller systems. In 2000, clusters containing 
about 10^ magnetic moments could be studied. This achievement rises the question 
whether further improvements might be possible. The fundamental limit of a SQUID 
is the quantum limit which correspond to a sensitivity of one magnetic moment 
for a SQUID with 1 ^m^. One might come close to this limit by using shunted 
SQUIDs [51]. Another possibility could be a reduction of the section of the micro- 
bridges [46]. Finally, the micro-SQUID technique could be improved by using 
superconducting materials with higher critical temperatures allowing measurements 
at higher temperatures. 



2f Mechanisms of Magnetization 
Reversal at Zero Kelvin 

As already briefly discussed in the introduction, for a sufficiently small magnetic 
sample it is energetically unfavorable to form a stable magnetic domain wall. The 
specimen then behaves as a single magnetic domain. For the smallest single domain 
particles \ the magnetization is expected to reverse by uniform rotation of magneti- 
zation (Sec. 2.1). For somewhat larger ones, non uniform reversal modes are more 
likely, for example the curling reversal mode (Sec. 2.2.1). For larger particles, mag- 
netization reversal occurs via a domain wall nucleation process starting in a rather 
small volume of the particle. For even larger particles, the nucleated domain wall can 
be stable for certain fields. The magnetization reversal happens then via nucleation 
and annihilation processes (Sec. 2.2.3). In these sections, we neglect temperature 
and quantum effects. 

The following section discusses in detail the uniform rotation mode which is used 
in many theories, in particular in Neel, Brown and Coffey's theory of magnetization 
reversal by thermal activation (Sect. 3) and in the theory of macroscopic quantum 
tunneling of magnetization (Sect. 4). 



2.1 MAGNETIZATION REVERSAL BY UNIFORM ROTATION 
(STONER-WOHLFARTH MODEL) 

The model of uniform rotation of magnetization, developed by Stoner and Wohl- 
farth [52J, and Neel [53], is the simplest classical model describing magnetization 
reversal. One considers a particle of an ideal magnetic material where exchange 
energy holds all spins tightly paraUel to each other, and the magnetization magnitude 
does not depend on space. In this case, the exchange energy is constant, and plays 
no role in the energy minimization. Consequently, there is competition only between 
the anisotropy energy of the particle and the effect of the appUed field. The original 
study by Stoner and Wohlfarth assumed only uniaxial shape anisotropy which is the 
anisotropy of the magnetostatic energy of the sample induced by its non-spherical 
shape. Thiaville has generalized this model for an arbitrary effective anisotropy 
which includes any magnetocrystalUne anisotropy and even surface anisotropy [54]. 



'in the theory of micromagnetism, the single domain state describes the state where the magnetization is 
perfectly aligned [2] whereas experimentalists mean often a state without domain wall. 
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In the simplest case of uniaxial anisotropy, the energy of a Stoner-Wohlfarth particle 
is given by 

E = KV sin^ (j) - hqMsVH cos(<?i - 6) , (2.1) 

where KV is the uniaxial anisotropy energy which depends on the shape of the 
particle, V is the volume of the particle, Ms is the spontaneous magnetization, H 
the magnitude of the appUed field, and cp and 9 are the angles of magnetization 
and applied field respectively, with respect to the easy axis of magnetization. The 
potential energy of (2.1) has two minima separated by an energy barrier. For given 
values of 9 and H, the magnetization Ues at an angle (p which locally minimizes 
the energy. This position can by found by equating to zero the first derivative with 
respect to </> of (2.1): dE/dcp = 0. The second derivative provides the criteria for 
maxima and minima. The magnetization reversal is defined by the minimal field 
value at which the energy barrier between the metastable minimum and the stable 
one vanishes, i.e., at dE/dcp = d'^E/dcjP' = 0. A short analysis yields the angular 
dependence of this field, called the switching field H^^, or in dimensionless units: 

77° 1 

(sin2/3 + cos2/3^) 

where = 2K/ (jioMs) is the anisotropy field. The angular dependence of h!^^ is 
plotted in Fig. 2.1. 

Contrary to h^^, the hysteresis loops cannot be expressed analytically and have to 
be calculated numerically. The result is seen in Fig. 2.2 showing the component 
of magnetization projected along the direction of the apphed field, i.e., Mh = 
Ms cos{(p — 6). Such loops are often called Stoner-Wohlfarth hysteresis loops. ^ 
The main advantage of this classical theory is that it is sufficiently simple to add some 
extra features to it, as presented in the following. 



2,1.1 Generalization of tlie Stoner-Wohlfarth model 

The original model of Stoner and Wohlfarth assumed only uniaxial shape anisotropy 
with one anisotropy constant (one second order term, see Eq. 2.1). This is sufficient 
to describe highly symmetric cases like a prolate spheroid of revolution or an infinite 
cyhnder. However, real systems are often quite complex and the anisotropy is a sum 
of mainly shape (magnetostatic), magnetocrystalline, magnetoelastic and surface 
anisotropy. One additional complication arises because the different contributions of 
anisotropics are often ahgned in an arbitrary way one with respect to each other. All 



^It is important to note that single particle measurement techniques do not measure this component Mh ■ 
For example for the micro-SQUlD technique, with the easy axis of magnetization in the plane of the 
SQUID and perpendicular to the current direction in the SQUID wire (Fig. 1.5), one measures a magnetic 
flux which is proportional to Ms sin (Fig. 2.3). 
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270° 

Fig. 2.1 Angular dependence of the S toner- Wohlfarth switching field = H^^/ H^, 
(Eq. 2.2). This curve is often called the 'Stoner-Wohlfarth astroid'. The cases (1) to (3) 
correspond to Eqs. 2.5 to 2.9 concerning the field dependence of the anisotropy barrier height. 




-1.5 -1 -0.5 0.5 1 1.5 
H/H, 



Fig. 2.2 Hysteresis loops of a S toner- Wohlfarth particle for different field angles 6. The 
component of magnetization along the direction of the applied field is plotted, i.e., Mn = 

Mscos((/« - 9). 
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these facts motivated a generalization of the Stoner-Wohlfarth model for an arbitrary 
effective anisotropy which was done by Thiaville [54, 55]. 

Similar to the Stoner-Wohlfarth model, one supposes that the exchange interaction 
in the cluster couples all the spins strongly together to form a giant spin whose 
direction is described by the unit vector m. The only degrees of freedom of the 
particle's magnetization are the two angles of orientation of m. The reversal of the 
magnetization is described by the potential energy: 

E{m, H) = Eo{m) - fioVM^m.H (2.3) 

where V and are the magnetic volume and the saturation magnetization of the 
particle respectively, H is the external magnetic field, and So(m) the magnetic 
anisotropy energy which is given by: 

Eo(m) = Esh3.pe{m) + EMcim) + EsnAa.ce{m) + EME{m) (2.4) 

-Eshapc is the magnetostatic energy related to the cluster shape. -Emc is the magne- 
tocrystaUine anisotropy (MC) arising from the coupling of the magnetization with 
the crystalline lattice, similar as in bulk. -Esurface is due to the symmetry breaking 
and surface strains. In addition, if the particle experiences an external stress, the 
volumic relaxation inside the particle induces a magneto-elastic (ME) anisotropy 
energy Eme- AU these anisotropy energies can be developed in a power series of 
m°m^m^ with p = a + b + c= 2,4,6,... giving the order of the anisotropy term. 
Shape anisotropy can be written as a biaxial anisotropy with two second order terms. 
Magnetocrystalline anisotropy is in most cases either unieixial (hexagonal systems) 
or cubic, yielding mainly second and forth order terms. Finally, in the simplest case, 
surface and magneto-elastic anisotropics are of second order. 

Thiaville proposed a geometrical method to calculate the particle's energy and to 
determine the switching field for all angles of the applied magnetic field yielding 
the critical surface of switching fields which is analogous to the Stoner-Wohlfarth 
astroid (Fig. 2.1). 

The main interest of Thiaville' s calculation is that measuring the critical surface of 
the switching field allows one to find the effective anisotropy of the nanoparticle. The 
knowledge of the latter is important for temperature dependent studies (Sect. 3) and 
quantum tunneling investigations (Sect. 4). Knowing precisely the particle's shape 
and the crystallographic axis allows one to determine the different contributions to 
the effective anisotropy. 

Thiaville' s calculation predicts also the field dependence of the energy barrier height 
AE close to switching (e = (l — H/H^?^) << 1) which is important to know for 
temperature dependent studies (Sects. 3 and 4). Three cases have to be distinguished: 

1 . In the majority of cases except the two following cases (see case 1 in Fig. 2.1): 



AE = 4KV^'-^s'/' = Eos^/' 
3 a/P 



(2.5) 
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where KV is the anisotropy energy constant, 7 is the angle of incidence 
between the local normal to the critical surface and the field sweeping direction, 
p is the radius of curvature of the focal curve [55] at H^^. It is important to 
emphasize that all these variables can be found experimentally by measuring 
the critical surface of the switching field. For uniaxial anisotropy, i.e. the 2D 
Stoner-Wohlfarth case, equation 2.5 becomes: 

/2\^/' |cos^|V3 

where 9 is the angle of the appUed field with respect to the easy axis of 
magnetization (Eq. 2.1). 



2. At glancing incidence (see case 2 in Fig. 2.1) with respect to the critical surface 
(7 = 7r/2), the power of e is different, yielding: 

l^E = E'gs^ (2.7) 

where Eq has been calculated only in the two dimensional case [54]. 

3. At a cusp point where 7 = 7r/2 (see case 3 in Fig. 2.1): 

AE = E%e'^ (2.8) 

where E"o has been calculated only in the two dimensional case [54]. In the 
case of uniaxial anisotropy, i.e. the 2D Stoner-Wohlfarth case), equation 2.8 
becomes: ^ 

AE = KVe'^ = Kv(^- (2.9) 

where iJa = -j^^ (Eq- 2.1). This equation is only valid for = and 7r/2, 
and it is valid for < H < Ha.. This famous result of Neel [9, 10] has often 
been wrongly used for assemblies of nanoparticles where it is very difficult to 
achieve the conditions — Q and 7r/2 [56]. Up to now, only the power 3/2 and 
2 (cases (1) and (3)) has been found by single particle measurements [36, 37] 
(see also Figs. 3.4 and 3.5). 



Another simple analytical approximation for the field dependence of the energy 
barrier AE{H) was derived numerically by Pfeiffer [57, 58]: 

AE^KV{l~H/H!!^y = Eoe% (2.10) 

where a = 0.86 + 1.14/1°^, and h'^^ is given by Eq. 2.2. This approximation is good 
for the intermediate field regime, i.e. for fields H not too close to H^, and not too 
small. 
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Fig. 2.3 Hysteresis loops of a nanocrystalline elliptic Co particle of 70 x 50 x 25 nm^. 
The dashed line is the prediction of the Stoner-Wohlfarth model of uniform rotation of 
magnetization. The deviations are due to non uniform magnetization states. 
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2.1.2 Experimental evidence for magnetization reversal by uniform rotation 

In order to demonstrate experimentally the uniform rotation mode, the angular de- 
pendence of the magnetization reversal has often been studied (see references in [2]). 
However, a comparison of theory with experiment is difficult because magnetic parti- 
cles often have a non uniform magnetization state which is due to rather complicated 
shapes and surfaces, crystalline defects, and surface anisotropy. In general, for many 
particle shapes the demagnetization fields inside the particles are non uniform lead- 
ing to non uniform magnetization states [2]. An example is presented in Fig. 2.3 
which compares typical hysteresis loop measurements of an elliptical Co particle, 
fabricated by electron beam lithography, with the prediction of the Stoner-Wohlfarth 
model (Fig. 2.3). Before magnetization reversal, the magnetization decreases more 
strongly than predicted because the magnetic configuration is not coUinear as in the 
Stoner-Wohlfarth model, but instead presents deviations mainly near the particle 
surface. The angular dependence of the switching field agrees with the Stoner- 
Wohlfarth model only for angles 9 < 30° where non-linearities and defects play a 
less important role [22, 59]. 

Studies of magnetization reversal processes in ultrathin magnetic dots with in-plane 
uniaxial anisotropy showed also switching fields which are very close to the Stoner- 
Wohlfarth model although magnetic relaxation experiments clearly showed that nu- 
cleation volumes are by far smaller than an individual dot volume [49]. These 
studies show clearly that switching field measurements as a function of the angles of 
the applied field cannot be taken unambiguously as a proof of a Stoner-Wohlfarth 
reversal. 

The first clear demonstration of the uniform reversal mode has been found with Co 
nanoparticles [36], and BaFeO nanoparticles [37], the latter having a dominant uniax- 
ial magnetocrystalUne anisotropy. The three dimensional angular dependence of the 
switching field measured on BaFeO particles of about 20 nm could be explained with 
the Stoner-Wohlfarth model taking into account the shape anisotropy and hexagonal 
crystalline anisotropy of BaFeO [38]. This explication is supported by temperature 
and time dependent measurements yielding activation volumes which are very close 
to the particle volume (Sect. 3). 

We present here the first measurements on individual cobalt clusters of 3 nm in 
diameter containing about a thousand atoms (Fig. 2.4) [39]. In order to achieve the 
needed sensitivity, Co clusters preformed in the gas phase are directly embedded in 
a co-deposited thin Nb film which is subsequently used to pattern micro-SQUIDs. 
A laser vaporization and inert gas condensation source is used to produce an intense 
supersonic beam of nanosized Co clusters which can be deposited in various matrices 
under Ultra High Vacuum (UHV) conditions. Due to the low energy deposition 
regime, clusters do not fragment upon impact on the substrate [60]. The niobium 
matrix is simultaneously deposited from a UHV electron gun evaporator leading to 
continuous films with a low concentration of embedded Co clusters [61]. These films 
are used to pattern planar micro-bridge-DC-SQUIDs by electron beam Uthography. 
The later ones allow us to detect the magnetization reversal of a single Co cluster for 
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Fig. 2.4 High Resolution Transmission Fig. 2.5 Scheme of a typical cluster shape 
Electron Microscopy observation along a with light-gray atoms belonging to the 1289- 
[110] direction of a 3 nm cobalt cluster ex- atoms truncated octahedron basis and dark- 
hibiting a/c.c-structure. gray atoms belonging to the (111) and (001) 

added facets. 

an applied magnetic field in any direction and in the temperature range between 0.03 
and 30 K (Sect. 1.2). However, the desired sensitivity is only achieved for Co clusters 
embedded into the micro-bridges where the magnetic flux coupling is high enough. 
Due to the low concentration of embedded Co clusters, we have a maximum of 5 
non-interacting particles in a micro-bridge which is 300 nm long and 50 nm wide. We 
can separately detect the magnetization switching for each cluster. Indeed they are 
clearly different in intensity and orientation because of the random distribution of the 
easy magnetization directions. The cold mode method (Sect. 1.2.5) in combination 
with the blind method (Sect. 1.2.6) allows us to detect separately the magnetic signal 
for each cluster. 

High Resolution Transmission Electron Microscopy observations showed that the 
Co clusters are well crystallized in a /c.c-structure (Fig. 2.4) with a sharp size 
distribution [61]. They mainly form truncated octahedrons (Fig. 2.5) [39]. 
Fig. 2.6 displays a typical measurement of switching fields in three dimensions 
of a 3 nm Co cluster at T = 35 mK. This surface is a three-dimensional picture 
directly related to the anisotropy involved in the magnetization reversal of the particle 
(Sect. 2.1). It can be reasonably fitted with the generalized Stoner and Wohlfarth 
model [54] (Sect. 2.1.1). We obtain the following anisotropy energy: 

Eo{rn)/v ~ —Kiml + K2m^ — i^4(m^,r7i^, + m'^,ml, + my,ml,) (2.11) 

where Ki and K2 are the anisotropy constants along z and x, the easy and hard 
magnetization axis, respectively. K4 is the forth order anisotropy constant and the 
{x'y'z') coordinate system is deduced from (xyz) by a 45° rotation around the z 
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axis. We obtained = 2.2 x 10^ J/m^, K2 = 0.9 x 10^ J/m^, and K4, = 0.1 x 10^ 
J/m^. The corresponding theoretical surface is showed in Fig. 2.7. Furthermore, we 
measured the temperature dependence of the switching field distribution (Sect. 3.3.2). 
We deduced the blocking temperature of the particle Tb « 14 K, and the number 
of magnetic atoms in this particle: N w 1500 atoms (Sect. 3.3.2). Detailed mea- 
surements on about 20 different particles showed similar three dimensional switch- 
ing field distributions with comparable anisotropy {Ki = (2.0 ± 0.3) x 10^ ilvc?, 
K2 = (0.8 ± 0.3) x 10^ J/m^, and = (0.1 ± 0.05) x 10^ J/m^) and size 
(N= 1500 ±200 atoms). 

In the following, we analyze various contributions to the anisotropy energy of the Co 
clusters. Fine structural studies using EXAFS measurements [6 1 ] were performed on 
500 nm thick niobium films containing a very low concentration of cobalt clusters. 
They showed that niobium atoms penetrate the cluster surface to almost two atomic 
monolayers because cobalt and niobium are miscible elements. Further magnetic 
measurements [61] on the same samples showed that these two atomic monolayers 
are magnetically dead. For this reason, we estimated the shape anisotropy of the 
typical nearly spherical deposited cluster in Fig. 2.5 after removing two atomic 
monolayers from the surface. By calculating all the dipolar interactions inside the 
particle assuming a bulk magnetic moment of = 1-7/iB. we estimated the shape 
anisotropy constants: Ki w 0.3 x 10^ J/m^ along the easy magnetization axis and 
K2 PS 0.1 X 10'' J/m^ along the hard magnetization axis. These values are much 
smaller than the measured ones which means that i^shape is not the main cause of the 
second order anisotropy in the cluster. 

The forth order term ii'4 = 0.1 x 10^ J/m^ should arise from the cubic magnetocrys- 
talline anisotropy in the f.c.c cobalt clusters. However, this value is smaller than the 
values reported in previous works [62, 63]. This might by due to the different atomic 
environment of the surface atoms in respect to that of hulk f.c.c. Co. Taking the 
value of the bulk [62, 63]: i^Tbuik = 1-2 x 10^ J/rn^ only for the core atoms in the 
cluster, we find i^MC ~ 0.2 x 10^ which is in reasonable agreement with our 
measurements. 

We expect that the contribution of the magnetoelastic anisotropy energy Kuc com- 
ing from the matrix-induced stress on the particle is also small. Indeed, using the 
co-deposition technique, niobium atoms cover uniformly the cobalt cluster creating 
an isotropic distribution of stresses. In addition, they can relax preferably inside the 
matrix and not in the particle volume since niobium is less rigid than cobalt. We 
believe therefore that only interface anisotropy -fl'surface can account for the exper- 
imentally observed second order anisotropy terms. Niobium atoms at the cluster 
surface might enhance this interface anisotropy through surface strains and magne- 
toelastic coupling. This emphasizes the dominant role of the surface in nanosized 
systems. 

In conclusion, the three dimensional switching field measurements of individual 
clusters giving access to its magnetic anisotropy energy. A quantitative understanding 
of the latter is still difficult but it seems that the cluster-matrix interface provides the 
main contribution to the magnetic anisotropy. Such interfacial effects could be 
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promising to control the magnetic anisotropy in small particles in order to increase 
their blocking temperature up to the required range for appUcations. 

2.1.3 Uniform Rotation with cubic anisotropy 

We have seen in the previous section that the magnetic anisotropy is often dominated 
by second order anisotropy terms. However for nearly symmetric shapes, forth order 
terms ^ can be comparable or even dominant with respect to the second order terms. 
Therefore, it is interesting to discuss further the features of forth order terms. We 
restrict the discussion to the 2D problem [64, 54] (see [55] for 3D). 
The reversal of the magnetization is described by equation 2.3 which can be rewritten 
in 2D: 

E{e) = Eo{6) - fiovM,{H^ cos{e) + Hy sin(e)) (2.12) 

where v and Af, are the magnetic volume and the saturation magnetization of the 
particle respectively, 6 is the angle between the magnetization direction and x, 
and Hy are the components of the external magnetic field along x and y, and Eq{6) 
is the magnetic anisotropy energy. The conditions of critical fields (dE/dO = and 
dE^/dO^ = 0) yield a parametric form of the locus of switching fields: 

^ f ■ rn.dE rn-.d^E\ 
= — sin 0) — +cos(6')— -r 

As an example we study a system with uniaxial shape anisotropy and cubic anisotropy. 
The total magnetic anisotropy energy can be described by: 

Eo{e) = vKi sin2(6» + 60) + VK2 sm^iO) 003^(9) (2.15) 

where Ki and K2 are anisotropy constants {Ki could be a shape anisotropy and K2 
the cubic crystalline anisotropy of a /c.c.-crystal). is a constant which allows to 
turn one anisotropy contribution with respect to the other one. Fig. 2.8 displays an 
example of a critical curve which can easily be calculated from Eqs. 2.13 to 2.15. 
When comparing the standard Stoner-Wohlfarth astroid (Fig. 2.1) with Fig. 2.8, we 
can realize that the critical curve can cross itself several times. In this case, the 
switching field of magnetization depends on the path followed by the applied field. 
In order to understand this point, let us follow the energy potential (2.12 and 2.15) 
when sweeping the applied field as indicated in Fig. 2.9. When the field is in A, 
the energy E has two minima and the magnetization is in the metastable potential 
well. As the field increases, the metastable well becomes less and less stable. Let 
us compare two paths, one going along A ^ Bi ^ C ^ D ^ E, the other over B2 



'For example, the forth order terms of f.c.c. magnetocrystalhne anisotropy. 



(2.13) 
(2.14) 
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Fig. 2.6 Top view and side view of the experimental three dimensional angular dependence 
of the switching field of a 3 nm Co cluster at 35 mK. This surface is symmetrical with respect 
to the Hx-Hy-plmie and only the upper part (iiqHz > T) is shown. Continuous lines on 
the surface are contour lines on which fioHz is constant. 




Fig. 2.7 Top view and side view of the theoretical switching field surface considering second 
and forth order terms in the anisotropy energy. 
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Fig. 2.8 Angular dependence of the switching field obtained from Eqs. 2.13 to 2.15 with 
i^i > and = -2/3Ki. The field is normaUzed by the factor 2Ki/{^ioMs). 




Fig. 2.9 Enlargement of angular depen- 
dence of the switching field of Fig. 2.8. Two 
possible paths of the applied field are indi- 
cated: starting from point A and going over 
the point Bi leads to magnetization reversal 
in E whereas going over the point B2 leads to 
reversal in D. 



Fig. 2.10 Scheme of the potential energy 
near the metastable state for different applied 
fields as indicated in Fig. 2.9. The ball repre- 
sent the state of the magnetization and the ar- 
row locate the appearing or disappearing well. 
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Fig. 2.11 Angular dependence of the switching field of a 3 nm Co cluster showing a strong 
influence of crystalline anisotropy. 



instead of Bi. Fig. 2.10 shows E in the vicinity of the metastable well for different 
field values along the considered paths (the stable potential well is not presented). 
One can realize that the state of the magnetization in C depends on the path followed 
by the field: going over Bi leads to the magnetization state in the left metastable 
well (1), whereas going over B2 leads to the right metastable well (2). The latter 
path leads to magnetization switching in D, and the former one to a switching in E. 
Note that a small magnetization switch happens when reaching Bi or B2. Point I is 
a supercritical bifurcation. 

The first measurement of such a field path dependence of the switching field were 
performed on single-domain FeCu nanoparticles of about 15 nm with a cubic crys- 
talline anisotropy and a small arbitrarily oriented shape anisotropy [65]. Fig. 2.11 
presents switching field measurements of a 15 rmi cobalt nanoparticle showing clearly 
a contribution of cubic crystalline anisotropy and the field path dependence of the 
switching field. 



2.2 NON UNIFORM MAGNETIZATION REVERSAL 

We have seen in the previous sections that for extremely small particles, magnetization 
should reverse by uniform rotation. For somewhat larger single domain particles, 
non uniform reversal modes are more likely. The simplest one is the curUng which 
is discussed in the following section. 
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2.2.1 Magnetization reversal by curling 

The simplest non uniform reversal mode is the curling reversal mode [66, 2]. The 
critical parameter is the exchange length A = ^fA/Mg, delimiting the region of 
uniform rotation and curling {A is the exchange constant). Therefore in the case of 
the size R > A (i? is for example the radius of a cylinder), magnetization reversal 
via curling is more favorable. In the following, we review the analytical result of an 
elUpsoid of rotation, which can by applied approximately to most of the shapes of 
nanoparticles or nanowires [67]. 

The variation of the switching field with the angle 9 (defined between the applied 
field and the long axis of the elUpsoid) is given by [68]: 

^sw = — ; > (2- 16) 

a1 s\v? 9 + 0^ cos^ 9 

where a^.z = "iNx.z — k/S"^, N^.z are the demagnetization factors, S = R/\, 
and R is the minor semi-axis of the ellipsoid. The parameter fc is a monotonically 
decreasing function of the aspect ratio of the ellipsoid. This function is plotted in 
Fig. 1 of Ref. [69]. The smallest and highest value of k is that for an infinite cylinder 
(k = 1 .079) and a sphere (k = 1 .379), respectively. 
For a long ellipsoid of rotation, the demagnetization factors are given by: 

Nz= ,\ ^L =ln(n-l + Vn2-l) , iV. = ^^, (2.17) 

where n is the ratio of the length to the diameter. 
For an infinite cylinder Eq. 2.16 becomes [66, 70]: 

Hl = ^^^^Ml^^=, (2.18) 
2 ^Kf + (1 + 2ht) cos2 9 

where ht = —1.079/5^. Eq. 2.18 is a good approximation for a very long elUpsoid 
of rotation. It is plotted in Fig. 2.12 for several radii of an infinite cylinder 
The case of uniform rotation of magnetization was generaUzed by Thiaville to an 
arbitrary anisotropy energy function and to three dimensions (Sect. 2.1.1). For the 
curling mode, this generalization is not possible. However, approximated calculations 
were proposed [71, 67, 72] and micro-magnetic simulations were performed [73]. 



2.2.2 Experimental evidence for magnetization reversal by curling 

We report here the first studies of isolated nanoscale wires with diameters smaller than 
100 rmi, for which single-domain states could be expected [35, 74]. The cylindrical 
geometry, with its large shape anisotropy, is well suited for comparison with theory. 
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Fig. 2.12 Angular dependence of the switching field of an infinite cylinder for several reduced 
cylinder radii S. For S < 1, the switching field is given by the uniform rotation mode 
(Sect. 2.1). 



Ni wires were produced by filling electrochemically the pores of commercially avail- 
able nanoporous track-etched polycarbonate membranes of thicknesses of 10 /xm. 
The pore size was chosen in the range of 30 to 100 imi [75, 76]. In order to place 
one wire on the SQUID detector, we dissolved the membrane in chloroform and put 
a drop on a chip of some hundreds of SQUIDs. Magnetization measurements were 
performed on SQUIDs with a single isolated wire (Fig. 1.1). 

The angular dependence of the switching field of wires with 45 and 92 nm in diameter 
are shown in Fig. 2.13. These measurements are in quantitative agreement with the 
curling mode (Eq. 2.16). Nevertheless, dynamical measurements showed a nucleation 
volume which is much smaller than the wire volume [35, 74] (Sect. 3.3.3). Therefore, 
we believe that the magnetization reversal starts close to curling instability but the 
nucleation happens in a small fraction of the wire only, then rapidly propagating 
along the whole sample. This picture is also in good agreement with micromagnetic 
simulations [73] and the micromagnetic model of Braun [77]. ^ 
The angular dependence of the switching field of Ni wires with larger diameters (270 
- 450 nm) were measured at room temperature by Lederman et al. [78]. Their results 
could roughly be explained by the curling mode. 



'*Note that the curling model only predicts an instability. There has never been a claim to describe what 
happens afterwards [2]. 
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Fig. 2.13 Angular dependence of the switching field of two Ni wires with a diameter of 45 
nm and 92 run, i.e. S =1.4 and 2.4. The switching fields are normalized by 125 mT and 280 
mT respectively. 



2.2.3 Magnetization reversal by nucleation and annihilation of domain walls 

For magnetic particles which have at least two dimetisiotis much larger than the 
domain wall width, the magnetization reversal may occur via nucleation/propagation 
and annihilation of one or several domain walls happening at two or more applied 
fields. We focus here on a 30-nm-thick elliptic Co particle defined by electron beam 
Uthography and Uft-off techniques out of sputtered thin films [79]. The Co film has a 
nanocrystalline structure leading to a magnetically soft material with a coercive field 
value of 3 mT at 4 K. Therefore we neglected the magneto-crystalline anisotropy. 
The nanofabricated particles have an elUptic contour with in-plane dimensions of 300 
nm X 200 nm and a thickness of 30 nm. 

In order to study the domain structure of our particles, we measure the angular 
dependence of hysteresis loops. Fig. 2.14 shows a typical hysteresis loop of an 
individual Co particle. The magnetic field is applied in the plane of the particle. The 
hysteresis loop is mainly characterized by two magnetization jumps. Starting from a 
saturated state, the first jump can be associated with domain wall nucleation and the 
second jump with domain wall annihilation. During these jumps, the magnetization 
switches in less than 100 /xs (our time resolution, see Sect. 1.2.4). The reversible 
central region of the hysteresis loops is evidence for the motion of the domain wall 
through the particle. 
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Fig. 2.14 Hysteresis loops at 0.1 K of the elliptic Co particle seen in the inset which shows 
the electron micrograph of the wire of a micro-bridge SQUID with the Co particle (300 nm 
X 200 nm x 30 nm). The in-plane field is applied along the long axis of the particles. The 
domain wall structure in an elliptical particle is also presented schematically as proposed by 
van den Berg. Arrows indicate the spin direction. The two small magnetization jumps near 
~ ±5mT might be due to the reversal of the center vortex of the domain structure. 



The simplest domain structure, showing such a hysteresis loop, has been proposed by 
van den Berg [80] in zero field and calculated for fields smaller than the saturation field 
by Bryant and Suhl [81]. ^ This domain structure has been observed experimentally 
on low anisotropy circular thin film disks (100 /im in diameter) using high-resolution 
Kerr techniques [82] . In zero field, the particle has a vortex-like domain wall as shown 
in Fig. 2.14. When a magnetic field is applied, this domain wall is pushed to the 
border of the particle. For higher fields the domain wall is annihilated and the particle 
becomes single domain. The main conditions of the model of van den Berg model are 
(i) that the magnetic material is very soft and (ii) that the system is two dimensional. 
The first condition is satisfied by our particles as they are made of a randomly oriented 
nanocrystallized Co, being very soft. The second condition is not quite well satisfied, 
however MFM measurements confirmed this domain structure [83]. Furthermore, 
we obtained similar results for thinner particles (10 and 20 nm). More complicated 
domain structures as proposed by van den Berg [79] may be excluded by the fact that 
similar Co particles of length smaller than 200 nm are single domain [22, 59]. 
After studying the magnetization reversal of individual particles, the question arises 
of which manner the properties of a macroscopic sample ai^e based on one-paiticle 



'Note that these models are 2D and neglect the domain wall width. Therefore, they can give only a 
qualitative description. 
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Fig. 2.15 Hysteresis loops of the magnetic moment of the array of 1.8 x 10^ Co particles 
(300 nm X 200 nm X 30 nm). The in-plane field is applied along the long axis of the particles. 



properties. In order to answer this question, we fabricated a sample consisting of 1 .8 
X 10*^ nearly identical elliptic Co particles of about the same dimensions and material 
as the individual particle studied above. These particles are placed on a Si substrate 
with a spacing of 2 /im. Because of this large spacing, dipole interactions between 
particles are negligible. Fig. 2. 15 shows the hysteresis loop of the array of Co particles 
when the field is applied parallel to the long axis of the particle. This hysteresis loop 
shows the same characteristics as the hysteresis loop of one particle (Fig. 2.14), i.e., 
nucleation and annihilation of domain walls. Because of switching field distributions 
mainly due to surface defects and a slight distribution of particle sizes, shapes, etc., 
the domain wall nucleation and annihilation are no longer discontinuous although 
still irreversible. They take place along continuous curves with a width of about 10 
mT. 



Influence of temperature on the 
magnetization reversal 

The thermal fluctuations of the magnetic moment of a single-domain ferromagnetic 
particle and its decay towards thermal equilibrium were introduced by Neel [9, 10] 
and further developed by Bean and Livingston [84, 85], and Brown [86, 87, 88]. The 
simplest case is an assembly of independent particles having no magnetic anisotropy. 
In the absence of an applied magnetic field, the magnetic moments are randomly 
oriented. The situation is similar to paramagnetic atoms where the temperature 
dependence of the magnetic susceptibility follows a Curie behavior, and the field 
dependence of magnetization is described by a Brillouin function. The only difference 
is that the magnetic moments of the particles are much larger than those of the 
paramagnetic atoms. Therefore, the quantum mechanical Brillouin function can be 
replaced by the classical limit for larger magnetic moments, namely the Langevin 
function. This theory is called superparamagnetism. The situation changes however 
as soon as magnetic anisotropy is present which establishes one or more preferred 
orientations of the particle's magnetization (Sect. 2). In the following, we present 
an overview over the simplest model describing thermally activated magnetization 
reversal of single isolated nanoparticles which is called the Neel-Brown Model. After 
a brief review of the model (Sect. 3.1), we present experimental methods to study 
the thermally activated magnetization reversal (Sect. 3.2). Finally, we discuss some 
applications of the Neel-Brown Model (Sect. 3.3). 



3.1 NEEL-BROWN MODEL OF THERMALLY ACTIVATED 
MAGNETIZATION REVERSAL 

In Neel and Brown's model of thermally activated magnetization reversal, a single 
domain magnetic particle has two equivalent ground states of opposite magnetization 
separated by an energy barrier which is due to shape and crystalline anisotropy. The 
system can escape from one state to the other by thermal activation over the barrier. 
Just as in the Stoner-Wohlfarth model, they assumed uniform magnetization and 
uniaxial anisotropy in order to derive a single relaxation time. Neel supposed further 
that the energy barrier between the two equilibrium states is large in comparison to 
the thermal energy k^T which justified a discrete orientation approximation [9, 10]. 
Brown criticized Neel's model because the system is not exphcitly treated as a 
gyromagnetic one [86, 87, 88]. Brown considered the magnetization vector in a 
particle to wiggle around an energy minimum, then jump to the vicinity of the other 
minimum, then wiggle around there before jumping again. He supposed that the 
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orientation of the magnetic moment may be described by a Gilbert equation with 
a random field term which is assumed to be white noise. On the basis of these 
assumptions. Brown was able to derive a Fokker-Planck equation for the distribution 
of magnetization orientations. Brown did not solve his differential equation. Instead 
he tried some analytic approximations and an asymptotic expansion for the case of 
the field parallel or perpendicular to the easy axis of magnetization. More recently, 
Coffey et al. [89, 90] found by numerical methods an exact solution of Brown's 
differential equation for uniaxial anisotropy and an arbitrary applied field direction. 
They also derived an asymptotic general solution for the case of large energy barriers 
in comparison to the thermal energy ksT. This asymptotic solution is of particular 
interest for single particle measurements and is reviewed in the following. 
For a general asymmetric bistable energy potential E = E{m, H) (Eq. 2.3) with the 
orientation of magnetization m = Ad /Ms (Ms is the spontaneous magnetization), 
H is the applied field, and with minima at ni and n2 separated by a potential 
barrier containing a saddle point at no (with the Hi coplanar), and in the case of 
l3{Eo - E,) » 1 where /? = l/fceT, and E, = E{n,, H), Coffey et al. showed 
that the longest relaxation time ' is given by the following equation which is valid in 
the intermediate to high damping limit (IHD) defined by a(3{Eo — Ei) > 1 [91]: 



-1 ^0 
T = 



27ra;o 

where wq and flo are the saddle and damped saddle angular frequencies: 



(3.1) 



^ J (0) (0) i\ 

^0 = TTV-Ci cy , (3.2) 



Ms 1 + a2 

wi and W2 are the well angular frequencies: 



1°) - + ^(4°)-ci°))2-4a-2ci°)4°) 



(3.3) 



'-. = ^V4^'c«, (3.4) 

withi = 1 and 2. c^^ andc2''^ (j = 0, 1, 2) are the coefficients in the truncated Taylor 
series of the potential at well and saddle points, i.e. the curvatures of the potential at 
well and saddle points. 7 is the gyromagnetic ratio, a = v^Ms is the dimensionless 
damping factor and v is the friction in Gilbert's equation (ohmic damping). 



The inverse of the longest relaxation time is determined by the smallest non- vanishing eigenvalue of the 
appropriate Fokker-Planck equation [89, 90]. All other eigenvalues can be neglected in the considered 
asymptotic limit of l3(Eo — Ei) » 1. 
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Whereas in the low damping limit (LD), defined by af3{Eo — Ei) < 1, the longest 
relaxation time is given by [92, 93]: 



-1 



a 



■l3(Eo-E2] 



(3.5) 



In this case, the energy dissipated in one cycle of motion in the well is very small in 
comparison to the thermal energy fce?"- 

Experimentally, relaxation is observed only if t is of the order of magnitude of 
the measuring time of the experiment. This implies for all known single particle 
measurement techniques that 0{Eo — Ei) >> 1, i.e. the asymptotic solutions (3.1) 
and (3.5) are always a very good approximation to the exact solution of Brown's 
Fokker-Planck equation [94]. Due to an applied field, (}{Eq - Ex) » (3{Eq - E2) 
(taking E2 as the metastable minimum) might be true. Then the first exponential in 
equations (3.1) and (3.5) can be neglected. 

Concerning the possible values of a, we remark that Uttle information is available. 

Typical values should be between 0.01 and 5 [8], meaning that in practice af3{Eo~Ei) 
can be >> 1, << 1, or 1. Thus the distinction between Eqs. (3.1) and (3.5) 
becomes important. 

Finally, we note that cY and (j = 0, 1,2) can be found experimentally by 
measuring the critical surface of the switching field and applying the calculation of 
Thiaville (Sect. 2.1.1) [55]. 



3.2 EXPERIMENTAL METHODS FOR THE STUDY OF THE 
NEEL-BROWN MODEL 

As discussed in the previous section, in the Neel-Brown model of thermally activated 
magnetization reversal, a single domain magnetic particle has two equivalent ground 
states of opposite magnetization separated by an energy barrier due to for instant 
shape and crystalline anisotropy. The system can escape from one state to the other 
either by thermal activation over the barrier at high temperatures, or by quantum 
tunneling at low temperatures (Sect. 4). At sufficiently low temperatures and at zero 
field, the energy barrier between the two states of opposite magnetization is much too 
high to observe an escape process. However, the barrier can be lowered by applying a 
magnetic field in the opposite direction to that of the particle's magnetization. When 
the applied field is close enough to the switching field at zero temperature -ffg^, 
thermal fluctuations are sufficient to allow the system to overcome the barrier, and 
the magnetization is reversed. 

In the following, we discuss three different experimental methods for studying this 
stochastic escape process which are called waiting time, switching field and telegraph 
noise measurements. 



38 



INFLUENCE OF TEMPERATURE ON THE MAGNETIZATION REVERSAL 




b) Switching field measurement 





Fig. 3.1 Schema of three methods for studying the escape from a potential well: waiting time 
and telegraph noise measurements give direct access to the switching time probability P{t), 
whereas switching field measurements yield histograms of switching fields 
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3.2.1 Waiting time measurements 

The waiting time method consists in measuring the probabiHty that the magnetization 
has not switched after a certain time. In the case of an assembly of identical and 
isolated particles, it corresponds to measurements of the relaxation of magnetization. 
However, in most particle assemblies, broad distributions of switching fields lead to 
logarithmic decay of magnetization, and the switching probability is hidden behind 
the unknown distributions functions [8]. For individual particle studies, waiting 
time measurements give direct access to the switching probability (Fig. 3.1). At 
a given temperature, the magnetic field H is increased to a waiting field iJ^ near 
the switching field H^^. Next, the elapsed time until the magnetization switches 
is measured. This process is repeated several hundred times, yielding a waiting 
time histogram. The integral of this histogram and proper normalization yields the 
probability that the magnetization has not switched after a time t. This probability is 
measured at different waiting fields i/w and temperatures in order to explore several 
barrier heights and thermal activation energies. 

According to the Neel-Brown model, the probability that the magnetization has not 
switched after a time t is given by: 

P{t) = e-*/- , (3.6) 

and T (inverse of the switching rate) can be expressed by an Arrhenius law of the 
form: 

T-i(e) Be°+''-ie-^^° , (3.7) 

where e = (l — H/H^^) and A, B, a and b depends on damping, temperature, 
energy barrier height (Eqs. 2.5 to 2.10, and 3.1 to 3.5), curvatures at well and saddle 
points, and reversal mechanism (thermal or quantum) (cf. Tab. 1 of Ref. [95]). For 
simplicity, experimentalists have often supposed a constant pre-exponential factor 
Tg-i instead of 

The adjustment of equation (3.6) to the measured switching probabilities yields a set 
of mean waiting times t^^ {H^,T). In order to adjust the Neel-Brown model to 
this set of data, we propose the following relation which can be found by inserting 
£= (l-ifw/^fsw)intoEq. (3.7): 



= H^^ 1 - 



\ In (tBs^+'-^) 



l/a\ 

(3.8) 



When plotting the values as a function of [Tin (rBe"+*-i)] all points 
should gather on a straight line (master curve) by choosing the proper value for the 
constants B, a and b (a and b should be given by Eqs. 2.5 to 2.10, and 3.1 to 3.5). A 
can be obtained from the slope of the master. 

The number of exploitable decades for t values is hmited for waiting time mea- 
surements: short-time (milliseconds) experiments are limited by the inductance of 
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the field coils ^ and long-time (minutes) studies by the stability of the experimental 
setup. Furthermore, the total acquisition time for a set of t~^{H^,T) is rather long 
(weeks). Thus a more convenient method is needed for single particle measurements 
which is the switching field method. 



3.2.2 Switching field measurements 

For single particle studies, it is often more convenient to study magnetization reversal 
by ramping the applied field at a given rate and measuring the field value as soon as 
the particle magnetization switches. Next, the field ramp is reversed and the process 
repeated. After several hundred cycles, switching field histograms are established, 
yielding the mean switching field (i?sw) and the width o-g^ (rms. deviation). Both 
mean values are measured as a function of the field sweeping rate and temperature 
(Fig. 3.1). 

From the point of view of thermally activated magnetization reversal, switching field 
measurements are equivalent to waiting time measurements as the time scale for the 
sweeping rate is typically more than 8 orders of magnitude greater than the time scale 
of the exponential prefactor which is in general around 10"^*^ s. We can therefore 
apply the Neel-Brown model described above. The mathematical transformation 
from a switching time probability (Eqs. 3.6 - 3.7) to a switching field probability 
was first given by Kurkijorvi [96] for the critical current in SQUIDs. A more general 
calculation was evaluated by Garg [95]. In many cases, the mean switching field 
(ifsw) can be approximated by the first two terms of the development of Garg [95]: 



1 



In 



vaA^-^/°- 



l/aN 



(3.9) 



where the field sweeping rate is given by u = dH/dt. The width of the switching 
field distribution CTsw can be approximated by the first term of Garg's development: 



y/6a 



l/a 



HLB 



(l-a)/a 



(3.10) 



In the case of a constant pre-exponential factor Tq ^, the calculation of (ffsw) and 
(Tsw is more simple and is given by equations (4) and (5) in Ref. [36], respectively. 
Similar to the waiting time measurements, a scaling of the model to a set of 
{Hsw{T,v)) values can be done by plotting the {Hsyf{T,v)) values as a function 



of 



Tin 



( g°w 



1/a 



AH points should gather on a straight line by choosing the 
proper value for the constants (a and b should be given by Eqs. 2.5 to 2.10, and 3.1 
to 3.5). 



solution to this problem might be a superposition of a constant applied field and a small pulse field. 
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The entire switching field distribution P{H) can be calculated iteratively by the 
following equation [96]: 



P{H) = T-\H)v-^ 



H 



P{H')dH' 



(3.11) 



3.2.3 Telegraph noise measurements 

In order to study the superparamagnetic state ^ of a single particle, it is simply 
necessary to measure the particle's magnetization as a function of time. We call 
this telegraph noise measurement as stochastic fluctuations between two states are 
expected. According to the Neel-Brown model, the mean time r spent in one state 
of magnetization is given by an Arrhenius law of the form of Eq. 3.7. As r increases 
exponentially with decreasing temperature, it is very unhkely that an escape process 
will be observed at low temperature. However, applying a constant field in direction 
of a hard axis (hard plane) of magnetization reduces the height of the energy barrier 
(Fig. 3.1). When the energy barrier is sufficiently small, the particle's magnetization 
can fluctuate between two orientations which are close to a hard axis (hard plane) 
of magnetization. The time spent in each state follows an exponential switching 
probability law as given by Eqs. 3.6 and 3.7 with a = 2 (Eq. 2.8). 



3.3 EXPERIMENTAL EVIDENCE FOR THE NEEL-BROWN MODEL 

The Neel-Brown model is widely used in magnetism, particularly in order to describe 
the time dependence of the magnetization of collections of particles, thin films and 
bulk materials. However until recently, all the reported measurements, performed on 
individual particles, were not consistent with the Neel-Brown theory. This disagree- 
ment was attributed to the fact that real samples contain defects, ends and surfaces 
which could play an important, if not dominant, role in the physics of magnetization 
reversal. It was suggested that the dynamics of reversal occurs via a complex path in 
configuration space, and that a new theoretical approach is required to provide a cor- 
rect description of thermally activated magnetization reversal even in single-domain 
ferromagnetic particles [19, 59]. Similar conclusions were drawn from numerical 
simulations of the magnetization reversal [97,98,99, 100, 101]. 



'At zero applied field, a single domain magnetic particle has two equivalent ground states of opposite 
magnetization separated by an energy barrier. When the thermal energy k^T is sufficiently high, the total 
magnetic moment of the particle can fluctuate thermally, like a single spin in a paramagnetic material. Such 
magnetic behavior of an assembly of independent single domain particles is called superparamagnetism 
[8, 84, 85]. 

'^Note that for a slightly asymmetric energy potential, one switching probability can be so long that 
two-level fluctuation becomes practically unobservable. 
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A few years later, micro-SQUID measurements on individual Co nanoparticles 
showed for the first time a very good agreement with the Neel-Brown model by 
using waiting time, switching field and telegraph noise measurements [36, 37, 39]. It 
was also found that sample defects, especially sample oxidation, play a crucial role 
in the physics of magnetization reversal. 

In the following subsections, we review some typical results concerning nanoparticles 
(Sect. 3.3.1), clusters (Sect. 3.3.2) and wires (Sect. 3.3.3). In Sect. 3.3.4, we point 
out the main deviations from the Neel-Brown model which are due to defects. 

3.3.1 Application to nanoparticles 

One of the important predictions of the Neel-Brown model concerns the exponential 
not-switching probability P{t) (Eq. 3.6) which can be measured directly via wait- 
ing time measurements (Sect. 3.2.1): at a given temperature, the magnetic field is 
increased to a waiting field which is close to the switching field. Then, the 
elapsed time is measured until the magnetization switches. This process is repeated 
several hundred times, in order to obtain a waiting time histogram. The integral of 
this histogram gives the not-switching probability P{t) which is measured at several 
temperatures T and waiting fields H„. The inset of Fig. 3.2 displays typical mea- 
surements of P{t) performed on a Co nanoparticle. All measurements show that 
P{t) is given by an exponential function described by a single relaxation time r. 
The validity of Eqs. 2.5 and 3.7) is tested by plotting the waiting field i?w as a 
function of [T/n(r/ro)]^/'^. ^ If the Neel-Brown model applies, all points should 
collapse onto one straight Une (master curve) by choosing the proper values for 
To. Fig. 3.2 shows that the data set r(iJw, T) falls on a master curve provided that 
To « 3 X 10~^s. The slope and intercept yield the values Eq = 214000 K and 
= 143.05 mT. The energy barrier Eq can be approximately converted to a thermally 
"activated volume" by using V = i?o/(A*o^'^s^sw) ~ (25 nm)^ which is very close 
to the particle volume estimated by SEM. This agreement is another confirmation 
of a magnetization reversal by uniform rotation. The result of the waiting time 
measurements are confirmed by switching field and telegraph noise measurements 
[36, 37]. The field and temperature dependence of the exponential prefactor tq is 
taken into account in [90]. 

3.3.2 Application to Co clusters 

Fig. 3.3 presents the angular dependence of the switching field of a 3 nm Co cluster 
measured at different temperatures. At 0.03 K, the measurement is very close to the 
standard Stoner-Wohlfarth astroid (Fig. 2.1). For higher temperatures the switching 
field becomes smaller and smaller. It reaches the origin at about 14 K yielding 
the blocking temperature Tb = 14 K of the cluster magnetization. Tb is defined 



'a = 3/2 because the field was applied at about 20° from the easy axis of magnetization (Eq. 2.5). 
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i:Tia(:T/T<,]|j?/3 

Fig. 3.2 Scaling plot of the mean switching time t{H^,T) for several waiting fields 

and temperatures (0.1 s < t{H^,T) < 60 s) for a Co nanoparticle. The scaling yields 
To ~ 3 X 10~^s. Inset: examples of the probability of not-switching of magnetization as 
a function of time for different applied fields and at 0.5 K. Full lines are data fits with an 
exponential function: P{t) = e"*/'" 



as the temperature for which the waiting time At becomes equal to the relaxation 

time T of the particle's magnetization at _ff = 0. Tb can be used to estimate the 
total number Ntot of magnetic Co atoms in the cluster. Using an Arrhenius-Uke law 
(Eq. 3.7) which can be written as: At = t = TQexp{KatNtot/k^T^), where t^^ is 
the attempt frequency typically between 10^" to 10^^ Hz [102], Kat is an effective 
anisotropy energy per atom and is the Boltzmann constant. Using the expression 
of the switching field at T = K and for ^ = 0: hqHsw = 2Ka,t/lJ.at = 0.3T 
(Fig. 3.3), the atomic moment jiat = ^-^ij-b, At = 0.01 s, tq = 10~^° s, and Tb = 
14 K, we deduce Ntot ~ 1500, which corresponds very well to a 3 nm Co cluster 
(Fig. 2.5). 

Fig. 3.4 presents a detailed measurement of the temperature dependence of the 
switching field at 0° and 45° and for three waiting times At. This measurement 
allows us to check the predictions of the field dependence of the barrier height. 
Eq. 3.9 predicts that the mean switching field should be proportional to T^/"^ where 
a depends on the direction of the appfied field (Eqs. 2.5 to 2.10): a = 2 for = 0° 
and 90°, and a = 3/2 for all other angles which are not too close to ^ = 0° and 90°. 
We found a good agreement with this model (Fig. 3.5). 
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Fig. 3.3 Temperature dependence of the switching field of a 3 nm Co cluster, measured in 
the plane defined by the easy and medium hard axes (Hy — plane in Fig. 2.6). The data 
were recorded using the blind mode method (Sect. 1.2.6) with a waiting time of the applied 
field of At = 0.1s. The scattering of the data is due to stochastic and in good agreement with 
Eq. 3.10. 
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Fig. 3.4 Temperature dependence of the switching field of a 3 nm Co cluster, measured at 
0° and 45°. The data were recorded using the blind mode method (Sect. 1.2.6) with different 
waiting time At of the applied field. The scattering of the data is due to stochastics and is in 
good agreement with Eq. 3.10. 




Fig. 3.5 Temperature dependence of the switching field of a 3 nm Co cluster as in Fig. 3.4 
but plotted as a function of (T/Tb)^/" with Tb = 14 K and a = 2 or 3/2 for 6* = 0° and 45°, 
respectively, and for three waiting times At. 
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3.3.3 Application to Ni wires 

Electrodepo sited wires (with diameters ranging from 40 to 100 nm and lengths up 
to 5000 nm, see Fig. 1.1) were studied [35, 74] using the micro-SQUID technique 
(Sect. 2.2.2). For diameter values under 50 nm, the switching probability as a 
function of time could be described by a single exponential function (Eq. 3.6). 
The mean waiting time r followed an Arrhenius law (Eq. 3.7) as proposed by the 
Neel-Brown model. Temperature and field sweeping rate dependence of the mean 
switching field could be described by the model of Kurkijarvi (Sect. 3.2.2) which is 
based on thermally assisted magnetization reversal over a simple potential barrier. 
This measurements allowed us to estimate an activation volume which was two orders 
of magnitude smaller than the wire volume. This confirmed the idea of the reversal 
of the magnetization caused by a nucleation of a reversed fraction of the cylinder, 
rapidly propagating along the whole sample. This result was also in good agreement 
with a micromagnetic model of Braun [77] . 

A pinning of the propagation of the magnetization reversal occurred for a few sam- 
ples, where several jumps were observed in the hysteresis curves. The pinning of 
a domain wall was probably due to structure defects. The dynamic reversal prop- 
erties of depinning were quite different from those of nucleation of a domain wall. 
For example, the probability of depinning as a function of time did not follow a 
single exponential law. A similar effect was also observed in single submicron Co 
particles having one domain wall [79], showing a domain wall annihilation process 
(Sect. 2.2.3). 

3.3.4 Deviations form the Neel-Brown model 

Anomalous magnetic properties of oxidized or ferrimagnetic nanoparticles have been 
reported previously by several authors [ 1 03, 1 04] . These properties are, e.g., the lack 
of saturation in high fields and shifted hysteresis loops after cooling in the presence 
of a magnetic field. These behaviors have been attributed to uncompensated surface 
spins of the particles and surface spin disorder [105, 106]. 

Concerning our single-particle studies, we systematically observed aging effects 
which we attribute to an oxidation of the surface of the sample, forming antifer- 
romagnetic CoO or NiO [59, 74]. We found that the antiferromagnetic coupling 
between the core of the particle or wire and its oxidized surface changed the dynamic 
reversal properties. For instant, we repeated the measurements of the magnetization 
reversal of a Ni wire two days after fabrication, six weeks after, and finally after 
three months [74]. Between these measurements, the wire stayed in a dry box. The 
quasi-static micro-SQUID measurements did reveal only small changes. The satu- 
ration magnetization measured after six weeks was unchanged and was reduced by 
one to two percent after three months. The angular dependence of the switching field 
changed also only slightly. The dynamic measurements showed the aging effects 
more clearly, as evidenced by 

• a non exponential probabiUty of not switching; 
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• an increase of the width of the switching field distributions; 

• a decrease of the activation energy. 

We measured a similar behavior on Uthographic fabricated Co piirticles with an 
oxidized border [59]. 

Fig. 3.6 presents the angular dependence of switching fields of a 3 nm Fe cluster 
having a slightly oxidized surface. Huge variation of the switching fields can be 
observed which might be due to exchange bias of frustrated spin configurations at 
the surface of the cluster (Fig. 3.7). 

We propose that the magnetization reversal of a ferromagnetic particle with an anti- 
ferromagnetic surface layer is mainly governed by two mechanisms which are both 
due to spin frustration at the interface between the ferromagnetic core and the antifer- 
romagnetic surface layer(s). The first mechanism may come from the spin frustration 
differing slightly from one cycle to another, thus producing a varying energy land- 
scape. These energy variations are less important at high temperatures when the 
thermal energy (A;bT) is much larger. However at lower temperature, the magne- 
tization reversal becomes sensitive to the energy variations. During the hysteresis 
loop the system chooses randomly a path through the energy landscape which leads 
to broad switching field distributions. A second mechanism may become dominant 
at high temperatures: the magnetization reversal may be governed by a relaxation of 
the spin frustration, hence by a relaxation of the energy barrier. This relaxation is 
thermally activated, i.e., slower at lower temperatures. 
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Fig. 3.6 Angular dependence of switching fields of a 3 nm Fe cluster having (probably) a 
slightly oxidized surface. Each point corresponds to on of the 10000 switching field measure- 
ment. The huge variations of the switching field might be due to exchange bias of frustrated spin 
configurations. However, quantum effects like those described in Sect. 4 are not completely 
excluded. 
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Fig. 3.7 Detail of the angular dependence of switching fields of a 3 nm Fe cluster having 
(probably) a slightly oxidized surface. Each point corresponds to one of the 3000 switching 
field measurement. Stochastic fluctuation between different switching field distributions are 
observed. The 'mean' hysteresis loop is shifted to negative fields. 



Magnetization reversal by 
quantum tunneling 

Studying the boundary between classical and quantum physics has become a very 
attractive field of research which is known as 'mesoscopic' physics (Fig. I.l). New 
and fascinating mesoscopic effects can occur when characteristic system dimensions 
are smaller than the length over which the quantum wave function of a physical quan- 
tity remains sensitive to phase changes. Quantum interference effects in mesoscopic 
systems have, until now, involved phase interference between paths of particles mov- 
ing in real space as in SQUIDs or mesoscopic rings. For magnetic systems, similar 
effects have been proposed for spins moving in spin space, such as magnetization 
tunneling out of a metastable potential well, or coherent tunneling between classically 
degenerate directions of magnetization [107, 108]. 

We have seen in the previous sections that the intrinsic quantum character of the 
magnetic moment can be neglected for nanoparticles with dimensions of the order of 
the domain wall width 5 and the exchange length A, i.e., particles with a collective 
spin of S* = 10^ or larger. However, recent measurements on molecular clusters with 
a collective spin of S = 10 suggest that quantum phenomena might be observed at 
larger system sizes with S » 1. Indeed, it has been predicted that macroscopic 
quantum tunneling of magnetization can be observed in magnetic systems with low 
dissipation. In this case, it is the tunneling of the magnetization vector of a single- 
domain particle through its anisotropy energy barrier or the turmeling of a domain 
wall through its pinning energy. These phenomena have been studied theoretically 
and experimentally [108]. 

The following sections review the most important results concerning the observed 
quantum phenomena in molecular clusters which are mesoscopic model systems to 
test quantum tunnehng theories and the effects of the environmental decoherence. 
Their understanding requires a knowledge of many physical phenomena, and are 
therefore particularly interesting for fundamental studies. We then focuses on Mag- 
netic Quantum Tunneling (MQT) studied in individual nanoparticles or nanowires. 
We concentrate on the necessary experimental conditions for MQT and review some 
experimental results which suggest that quantum effects might even be important in 
nanoparticles with S = 10^ or larger. 
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Fig. 4.1 Schematic view of the magnetic core of the Fes cluster. The oxygen atoms are 
black, the nitrogen gray, and carbon atoms are white. For the sake of clarity only the hydrogen 
atoms that are exchanged with deuterium are shown as small spheres (Sect. 4.2.4). The arrows 
represent the spin structure of the ground state S = 10 as experimentally determined through 
polarized neutron diffraction experiments [109]. The exact orientation of easy, medium and 
hard axis of magnetization (Fig. 4.2) can be found in [110]. 



4.1 QUANTUIVI TUNIVELING OF IMAGISfETIZATION EV IMOLECULAR 
CLUSTERS 

Magnetic molecular clusters are the final point in the series of smaller and smaller 
units from bulk matter to atoms. Up to now, they have been the most promising 
candidates for observing quantum phenomena since they have a well defined structure 
with well characterized spin ground state and magnetic anisotropy. These molecules 
can be regularly assembled in large crystals where all molecules often have the 
same orientation. Hence, macroscopic measurements can give direct access to single 
molecule properties. The most prominent examples are a dodecanuclear mixed- 
valence manganese-oxo cluster with acetate ligands, short Mni2 acetate [111], and 
an octanuclear iron(III) oxo-hydroxo cluster of formula [Fe802(OH)i2(tacn)6]^"'" 
where tacn is a macrocyclic ligand, short Feg (Fig. 4.1) [112]. Both systems have 
a spin ground state of 5 = 10, and an Ising-type magneto-crystalline anisotropy, 
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which stabilizes the spin states with m = ±10 and generates an energy barrier for 
the reversal of the magnetization of about 67 K for M1112 acetate and 25 K for Feg. 
Thermally-activated quantum tunneling of the magnetization has first been evidenced 
in both systems [113, 114, 115, 116]. Theoretical discussion of this assumes that 
thermal processes (principally phonons) promote the molecules up to high levels 
with small quantum numbers |m|, not far below the top of the energy barrier, and the 
molecules then tunnel inelastically to the other side. Thus the transition is almost 
entirely accomplished via thermal transitions and the characteristic relaxation time is 
strongly temperature dependent. An alternative explication was also presented [117]. 
For Fes however the relaxation time becomes temperature independent below 0.36 
K [1 16, 118] showing that a pure tunneling mechanism between the only populated 
ground states m = ±5* = ±10 is responsible for the relaxation of the magnetization. 
On the other hand in the Mni2 acetate system one sees temperature independent 
relaxation only for strong applied fields and below about 0.6 K [1 19, 120]. During 
the last years, several new molecular magnets were presented (see for instance [121, 
122, 123, 124]) which show also tunneling at low temperatures. 
The following subsections review the most appealing results concerning the Feg 
system with can be seen as an ideal 'model molecule' to study quantum phenomena 
in magnetic nanostructures. We stress that the tunneling in large spins is remarkable 
because it does not show up at lowest orders of perturbation theory. 
All measurements on Feg were performed using an array of micro-SQUIDs (Sect. 1 .2.7), 
The high sensitivity of this magnetometer allows us to study single Fes crystals [125] 
of sizes of the order of 10 to 500 /xm. For ac-susceptibility measurements and 
magnetization measurements at T > 6 K, we used a home built Hall probe magne- 
tometer [26, 126]. It works in the temperature range between 0.03 K and 100 K, for 
frequencies between 1 Hz and 100 kHz. 

After discussing the magnetic anisotropy of Fes, we present the observed quantum 
phenomena. The discussions of the following sections neglect environmental deco- 
herence effects for the sake of simplicity. In Sects. 4.2, we focus on effects of the 
environment (dipolar coupling, nuclear spins, and temperature) onto the turmeUng. 
This review should help to set up a complete theory which describes real magnetic 
quantum systems. 

4,1.1 Magnetic anisotropy in Fes 

The octanuclear iron(III) oxo-hydroxo cluster of formula [Fe802(OH)i2(tacn)6]*"'" 
where tacn is a macrocyclic ligand, short Feg (Fig. 4.1) was first synthesized by 
Wieghardt ef flZ. in 1984 [125]. Four central iron (III) ions with 5 = 5/2 are bridged 
by two 0x0 groups. The other four iron (III) ions are bridged by hydroxo groups 
to the central iron ions in an almost planar arrangement, as shown in Fig. 4.1. The 
clusters have approximate D2 symmetry but crystallize in the triclinic system [110]. 
Feg has a S = 10 ground state which originates from antiferromagnetic interactions 
that do not give complete compensation of the magnetic moment [127]. Spin- 
orbital moments can be neglected as the magnetic ions are in an 'orbital singlet' as 
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a result of Hund's rules. The spin structure of the ground state schematized by the 
arrows in Fig. 4.1 has been recently confirmed by a single-crystal polarized neutron 
investigation that provided a magnetization density map of the cluster [109]. 

The simplest model describing the spin system of Fcs molecular clusters (called the 
giant spin model) has the following Hamiltonian 1 112]: 

H = -DSl + E{Sl-Sl)+gi^Bl^oS-H . (4.1) 

Sx, Sy, and are the three components of the spin operator, D and E are the 
anisotropy constants which were determined via HF-EPR {D/k^ « 0.275 K and 
E/k-QK. 0.046 K [1 12]), and the last term of the Hamiltonian describes the Zeeman 
energy associated with an applied field H. This Hamiltonian defines a hard, medium, 
and easy axes of magnetization in x, y and z directions, respectively (Fig. 4.2). It has 
an energy level spectrum with (25 + 1) = 21 values which, to a first approximation, 
can be labeled by the quantum numbers m = —10, —9, 10 choosing the z-axis 
as quantization axis. The energy spectrum, shown in Fig. 4.3, can be obtained by 
using standard diagonalisation techniques of the [21 x 21] matrix describing the spin 
Hamiltonian S* = 10. At H = 0, the levels to = ±10 have the lowest energy. 
When a field is applied, the energy levels with m < — 2 increase, while those 
with TO > 2 decrease (Fig. 4.3). Therefore, energy levels of positive and negative 
quantum numbers cross at certain fields H^- It turns out that for Feg the levels cross 
at fields given by /Uo-ffz ~ nx 0.22 T, with n = 1,2, 3.... The inset of Fig. 4.3 
displays the details at a level crossing where transverse terms containing Sx or Sy 
spin operators turn the crossing into an 'avoided level crossing'. The spin S is 'in 
resonance' between two states when the local longitudinal field is close to an avoided 
level crossing. The energy gap, the so-called 'tunnel spitting' A, can be tuned by 
an applied field in the xy— plane (Fig. 4.2) via the SxH^ and SyHy Zeeman terms 
(Sect. 4.1.3). 

The effect of these avoided level crossings can be seen in hysteresis loop mea- 
surements (Fig. 4.4). When the applied field is near an avoided level crossing, the 
magnetization relaxes faster, yielding steps separated by plateaus. As the temperature 
is lowered, there is a decrease in the transition rate due to reduced thermal-assisted 
tunneling. A similar behavior was observed in Mni2 acetate clusters [113, 114, 115] 
where equally separated steps were observed at i?z « n x 0.45 T. The main difference 
between both clusters is that the hysteresis loops of Feg become temperature inde- 
pendent below 0.36 K whereas measurements on Mni2 acetate indicate a temperature 
independence only for strong applied fields and below 0.6K [119, 120]. 
Another important difference is that the step heights {i.e. the relaxation rates) change 
periodically when a constant transverse field is applied (Fig. 4.5). It is the purpose 
of the next subsections to present a detailed study of this behavior. 
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Fig. 4.2 Unit sphere showing degenerate minima A and B which are joined by two tunnel 
paths (heavy lines). The hard, medium, and easy axes are taken in x-, y- and 2;-direction, 
respectively. The constant transverse field , „„., for tunnel splitting measurements is applied 
in the a;y-plane at an azimuth angle (p. At zero applied field H = 0, the giant spin reversal 
results from the interference of two quantvmi spin paths of opposite direction in the easy 
anisotropy ?yz-planc. For transverse fields in direction of the hard axis, the two quantum spin 
paths are in a plane which is parallel to the j/z-plane, as indicated in the figure. By using Stokes 
theorem it has been shown [128] that the path integrals can be converted in an area integral, 
giving that destructive interference, that is a quench of the tunneling rate, occurs whenever the 
shaded area is kn/S, where k is an odd integer. The interference effects disappear quickly 
when the transverse field has a component in the ^/-direction because the tunneUng is then 
dominated by only one quantimi spin path. 
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Fig. 4.3 Zeeman diagram of the 21 levels of the = 10 manifold of Fes as a function of 
the field applied along the easy axis (Eq. 4.1). From bottom to top, the levels are labeled with 
quantum numbers m = ±10, ±9, ...0. The levels cross at fields given by /io^^z ~ nx 0.22 
T, with n = 1, 2, 3.... The inset displays the detail at a level crossing where the transverse 
terms (terms containing or/and Sy spin operators) turn the crossing into an avoided level 
crossing. The greater the tunnel splitting A, the higher the tunnel rate. 



4.1.2 Landau-Zener tunneling in Fcg 

The non-adiabatic transition between the two states in a two-level system has first 
been discussed by Landau, Zener and Stiickelberg [130, 131, 132]. The original 
work by Zener concentrates on the electronic states of a bi-atomic molecule, while 
Landau and Stiickelberg considered two atoms which undergo a scattering process. 
Their solution of the time dependent Schrddinger equation of a two-level system 
could be applied to many physical systems and it became an important tool for 
studying tunneling transitions. The Landau-Zener model has also been applied to 
spin tunneling in nanoparticles and clusters [133, 134, 135, 136, 137, 138]. The 
tunneling probability P when sweeping the longitudinal field at a constant rate 
over an avoided energy level crossing (Fig. 4.6) is given by: 



Here, m and m' are the quantum numbers of the avoided level crossing, dHz/dt is 
the constant field sweeping rates, 5 ~ 2, /Ub the Bohr magneton, and h is Planck's 
constant. 

With the Landau-Zener model in mind, we can now start to understand qualitatively 
the hysteresis loops (Fig. 4.4b). Let us start at a large negative magnetic field H^. 





(4.2) 



2Ugij,B\m — m'\fiQdHz/dt 
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Fig. 4.4 Hysteresis loops of a single crystal of Peg molecular clusters: (a) at different 
temperatures and a constant sweeping rate dHz/dt = 0.014 T/s, and (b) at 0.04 K and 
different field sweeping rates. The loops display a series of steps, separated by plateaus. As the 
temperature is lowered, there is a decrease in the transition rate due to reduced thermal assisted 
tunneling. The hysteresis loops become temperature independent below 0.35 K, demonstrating 
quantum tunneUng at the lowest energy levels. The resonance widths at small fields of 
about 0.05 T are mainly due to dipolar fields between the molecular clusters [118, 129]. 
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Fig. 4.5 Hysteresis loops measured along Hz in the presence of a constant transverse field 
at 0.04 K. Insets: enlargement around the field Hz = 0. Notice that the sweeping rate is ten 
times slower for the measurements in the insets than that of the main figures. 
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Fig. 4.6 Detail of the energy level diagram near an avoided level crossing, m, and m' are the 
quantum numbers of the energy level. Pm,m' is the Landau-Zener tunnel probability when 
sweeping the applied field from the left to the right over the anti-crossing. The greater the gap 
A and the slower the sweeping rate, the higher is the tunnel rate (Eq. 4.2). 



At very low temperature, all molecules are in the m = — 10 ground state. When the 
applied field Hz is ramped down to zero, all molecules will stay in the m = — 10 
ground state. When ramping the field over the A_io.io-region at ~ 0, there 
is a Landau-Zener tunnel probability -P_io,io to tunnel from the m = —10 to the 
m = 10 state. P_io,io depends on the sweeping rate (Eq. 4.2), i.e. the slower the 
sweeping rate, the larger _P_io.io- This is clearly demonstrated in the hysteresis loop 
measurements showing larger steps for slower sweeping rates (Fig. 4.4b). When the 
field Hz is now further increased, there is a remaining fraction of molecules in the 
m = — 10 state which became a metastable state. The next chance to escape from this 
state is when the field reaches the A_io,9-region. There is a Landau-Zener tunnel 
probability P-10,9 to tunnel from the m = —10 to the m — 9 state. As m = 9 is 
an excited state, the molecules in this state desexcite quickly to the m = 10 state by 
emitting a phonon. An analogous procedure happens when the applied field reaches 
the A_io.io-n-regions (n = 2, 3, . . . ) until all molecules are in the m = 10 ground 
state, i.e. the magnetization of all molecules is reversed. As phonon emission can 
only change the molecule state by Am = 1 or 2, there is a phonon cascade for higher 
applied fields. ^ 

In order to apply quantitatively the Landau-Zener formula (Eq. 4.2), we first saturated 
the crystal of Fes clusters in a field ofH^ = — 1 .4 T, yielding an initial magnetization 



Phonon induced transitions with | Am| > 2 are very small [139, 140, 141]. 



I m'> 
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Min = — Mg. ^ Then, we swept the applied field at a constant rate over one of the 
resonance transitions and measured the fraction of molecules which reversed their 
spin. This procedure yields the tunnehng rate P_io,io-ti and thus the tunnel spUtting 
A-io,io-n (Eq. 4.2) with n = 0, 1, 2, . . . . 

For very small turmeling probabilities P_io,io-raj we did multiple sweeps over the 
resonance transition. The magnetization M after N sweeps is given by: 

M{N) « Meq + (Min - Meq)e-'=^-^°.i°-''^ = M^q + (Min - Meq)e-" (4.3) 

Here Min is the initial magnetization, Mgq{Hz) is the equilibrium magnetization, 
TV = -J is the number ofsweeps over the level crossing, r = fcP_io,io-n^^^ 
is the overall Landau-Zener transition rate, k = 2 for n = and k = 1 for n = 1,2, 
. . . , and A is the amplitude of the ramp-field. ^ We have therefore a simple tool to 
obtain the tunnel splitting by measuring P_io,io-n, or M{N) for P_io,io-n << 1- 
We first checked the predicted Landau-Zener sweeping field dependence of the tun- 
neling rate. This can be done, for example, by plotting the relaxation of magnetization 
as a function oft = Njj^j^- The Landau Zener model predicts that all measure- 
ments should fall on one curve which was indeed the case for sweeping rates between 
1 and 0.001 T/s (Fig. 4.7) for the m = ±10 transition. The deviations at lower 
sweeping rates, are mainly due to the hole-digging mechanism [129] which slows 
down the relaxation (see Sect. 4.2.2). "* In the ideal case, we should find an exponen- 
tial curve (Eq. 4.3). However, we found clear deviations from the exponential curve 
(Fig. 4.8). which might be due to molecules with different amounts of nuclear spins. 
For example, two percent of natural iron has a nuclear spin, i.e. about 10 percent 
of Fes has at least one nuclear spin on the iron. This interpretation is supported by 
measurements on isotopically substituted Feg samples (Fig. 4.8). 
Another way of checking the Landau-Zener sweeping field dependence of the tun- 
nehng rate is presented in Fig. 4.9 showing a sweeping rate independent A_io,io 
between 1 and 0.001 T/s. ^ The measurements on isotopically substituted Fes sam- 
ples show a small dependence of A^,™' on the hyperfine coupling (Fig. 4.9). Such 
an effect has been predicted for a constant applied field by Tupitsyn et al. [142], and 
for a ramped field by Rose [136]. Further details are presented in Sect. 4.2.4. 
We also compared the turmeling rates found by the Landau-Zener method with those 
found using a square-root decay method which was proposed by Prokof'ev and 
Stamp [143], and found a good agreement [129, 144] (Sect. 4.2.1). 
Our measurements showed for the first time that the Landau-Zener method is partic- 
ularly adapted for molecular clusters because it works even in the presence of dipolar 



^In order to avoid heating problems for measurements of A forn > 1, we started in a thermally annealed 
sample with Min = 0.95Ms instead of Min = -A^s or Min = 0. 

'We supposed here that the forth and back sweeps give the same tunnel probability. This is a good 
approximation for P << 1 where next nearest neighbor (molecule) effects can be neglected. 
'^Roughly speaking, at very low field sweeping rates internal fields change faster than the external field. 
^Recent measurements confirmed the good agreement up to 30 T/s. 
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Fig. 4.7 Scaling plot for the Landau-Zener method showing the predicted field sweeping 
rate dependence for 1 T/s to 1 mT/s. Each point indicates the magnetization after a field sweep 
over the m = ±10 resonance. The dotted lines are guides for the eyes. For comparison, the 
figure displays also a relaxation curve at a constant field H = (Fig. 4.5) which shows much 
slower relaxation [129]. 



fields which spread the resonance transition provided that the field sweeping rate 
is not too small. Furthermore, our measurements show a small but clear influence 
of the hyperfine coupling which should be included in a generaUzed Landau-Zener 
model [136]. 

4.1.3 Oscillations of tunnel splitting 

An applied field in the plane can tune the tunnel sphttings Am,m' via the 
and Sy spin operators of the Zeeman terms which do not commute with the spin 
Hamiltonian. This effect can be demonstrated by using the Landau-Zener method 
(Sect. 4.1.2). Fig. 4.10 presents a detailed study of the tunnel sphtting A±io at the 
tunnel transition between m = ±10, as a function of transverse fields applied at 
different angles y?, defined as the azimuth angle between the anisotropy hard axis 
and the transverse field (Fig. 4.2). For small angles (p the tunneling rate oscillates 
with a period of ca. 0.4 T, whereas no oscillations showed up for large angles 
(f [47]. In the latter case, a much stronger increase of A±io with transverse field 
is observed. The transverse field dependence of the tunnehng rate for different 
resonance conditions between the state m = — 10 and (10 — n) can be observed by 
sweeping the longitudinal field around hqHz = n x 0.22 T with n = 0, 1, 2, ... The 
corresponding tunnel splittings A_io,io-n oscillate with almost the same period of 
ca. 0.4 T (Fig. 4.10). In addition, comparing quantum transitions between m = — 10 
and (10 — n), with n even or odd, revealed a parity (or syimnetry) effect which 
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Fig. 4.8 Scaling plot for the Landau-Zener method showing the predicted field sweeping 
rate dependence for 1 T/s to 1 mT/s similar to Fig. 4.7 but for three isotopically substituted 
Feg samples. Further details are presented in Sect. 4.2.4 




Fig. 4.9 Field sweeping rate dependence of the tunnel splitting A^io.io measured by a 
Landau-Zener method for three Fes samples, for Hx = 0. The Landau-Zener method works 
in the region of high sweeping rates where A_io,io is sweeping rate independent. Note that 
the differences of A_io.io between the three isotopically substituted samples are rather small 
in comparison to the oscillations in Fig. 4.10. 
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fjiagri^tj iroirisiij^rs^^ fiSfs! CD- 
Fig. 4.10 Measured tunnel splitting A as a function of transverse field for (a) several azimuth 

angles ip atm = ±10, and (b) f x 0°, and for quantum transition between m = — 10 and 
(10 — n). Note the parity effect which is analogous to the suppression of tunneling predicted 
for half integer spins. It should also be mentioned that internal dipolar and hyperflne fields 
hinder a quench of A which is predicted for an isolated spin (Figs. 4.11 and 4.14). 



is analogous to the Kramers' suppression of tunneling predicted for half integer 
spins [145, 146]. This behavior has been observed for n = to 4. * A similar 
strong dependence on the azimuth angle if was observed for all the resonances. 



*The tunneling rate were too fast for n > 4. 
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4.1.3.1 Semi-classical descriptions Before showing that the above results can be 
derived by an exact numerical calculation using the quantum operator formalism, it 
is useful to discuss semi-classical models. The original prediction of oscillation of 
the tunnel sphtting was done by using the path integral formalism [147]. Here [128], 
the oscillations are explained by constructive or destructive interference of quantum 
spin phases (Berry phases) of two tunnel paths (instanton trajectories) (Fig. 4.2). 
Since our experiments were reported, the Wentzel-Kramers-Brillouin theory has 
been used independently by Garg [148] and Villain and Fort [149]. The surprise is 
that although these models [128, 148, 149] are derived semiclassically, and should 
have higher order corrections in 1/S, they appear to be exact as written! This has 
first been noted in Refs. [148, 149] and then proven in Ref . [150] Some extensions or 
alternative explications of Garg's result can be found in Refs. [151, 152, 153, 154]. 
The period of oscillation is given by [128]: 

AH = —y/2E(E + D) (4.4) 

where D and E are defined in Eq. 4.1. We find a period of oscillation of AH = 0.26 T 
for D = 0.275 K and = 0.046 K as in Ref. [112]. This is somewhat smaller than 
the experimental value of ca. 0.4 T. We believe that this is due to higher order terms 
of the spin Hamiltonian which are neglected in Garg's calculation. These terms can 
easily be included in the operator formaUsm as show in next subsection. 

4.1.3.2 Exact numerical diagonalisation In order to quantitatively reproduce the 
observed periodicity we included forth order terms in the spin Hamiltonian (Eq. 4. 1) 
as recently employed in the simulation of inelastic neutron scattering measurements 
[155, 156] and performed a diagonalisation of the [21 x 21] matrix describing the 
S' = 10 system. For the calculation of the tunnel splitting we used D — 0.289 K, 
E = 0.055 K(Eq. 4.1) and the forth order terms as defined in [155] withB^ = 0.72x 
10^6 K, Bl = l.Olx 10"^ K, Bi ^ -0.43X 10"'' K, which are close to the values 
obtained by EPR measurements in [110] and neutron scattering measurements [156]. 
The calculated tunnel splittings for the states involved in the tunneling process at 
the resonances n = 0, 1, and 2 are reported in Fig. 4.11, showing the oscillations 
as well as the parity effect for odd resonances. The calculated tunneling splitting is 
however ca. 1.5 times smaller than the observed one. This small discrepancy could 
be reduced by introducing higher order terms. We believe that this is not relevant 
because the above model neglects for example the influence of nuclear spins which 
seems to increase the measured (effective) tunnel splittings (Fig. 4.9 and Sect. 4.2.4). 
Our choice of the forth order terms suppresses the oscillations of large transverse 
fields (Fig. 4.11). This region could not be studied in the current set-up. Future 
measurements should focus on the higher field region in order to find a better effective 
Hamiltonian. 
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Fig. 4.11 Calculated tunnel splitting A as a function of transverse field for (a) quantum 
transition between m = ±10 at several azimuth angles (f, and (b) quantum transition between 
m = —10 and (10 — n) at ip = 0° (Sect. 4.1.3.2). The forth order terms suppress the 
oscillations of A for large transverse fields | Hx \ . 
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4.2 ENVIRONMENTAL DECOHERENCE EFFECTS IN MOLECULAR 
CLUSTERS 

At temperatures below 0.36 K, Fcg molecular clusters display a clear crossover 
from thermally activated relaxation to a temperature independent quantum regime, 
with a pronounced resonance structure of the relaxation time as a function of the 
external field (Sect. 4.1.1). It was surprising however that the observed relaxation of 
the magnetization in the quantum regime was found to be non-exponential and the 
resonance width orders of magnitude too large [116, 118]. The key to understand this 
seemingly anomalous behavior involves the hyperfine fields as well as the evolving 
distribution of the weak dipole fields of the nanomagnets themselves L143J. Both 
effects showed to be the main source of decoherence at very low temperature. At 
higher temperatures, phonons are another source of decoherence. 
In the following sections, we focus on the low temperature and low field Umits, 
where phonon-mediated relaxation is astronomically long and can be neglected. In 
this limit, the m = ±5 spin states are coupled due to the tunneling splitting A±s 
which is about 10"^ K for Feg (Sect. 4.1.3) and 10"" K for Mni2 [157] with S 
= 10. In order to tunnel between these states, the longitudinal magnetic energy 
bias ^ = .9MB'S'-ffiocai due to the local magnetic field -ffiocai on a molecule must be 
smaller than A±s implying a local field smaller than 10~* T for Fes clusters. Since 
the typical intermolecular dipole fields are of the order of 0.05 T, it seems at first that 
almost all molecules should be blocked from tunneling by a very large energy bias. 
Prokof 'ev and Stamp have proposed a solution to this dilemma by proposing that fast 
dynamic nuclear fluctuations broaden the resonance, and the gradual adjustment of 
the dipole fields in the sample caused by the tunneling, brings other molecules into 
resonance and allows continuous relaxation [143]. Some interesting predictions are 
briefly reviewed in the following section. 

4.2.1 Prokof 'ev-Stamp theory 

Prokof ev and Stamp were the first who realized that there are localized couplings 
of environmental modes with mesoscopic systems which cannot be modeled with 
an 'oscillator bath' model [158] describing delocalized environmental modes such 
as electrons, phonons, photons, etc. They found that these localized modes such as 
nuclear and paramagnetic spins are often strong and described them with a spin bath 
model [159]. We do not review this theory ^ but focus on one particular application 
which is interesting for molecular clusters [143]. Prokof'ev and Stamp showed that 
at a given longitudinal appUed field H^, the magnetization of a crystal of molecular 
clusters should relax at short times with a square-root time dependence which is due 



For a review, see [160]. 
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to a gradual modification of the dipole fields in the sample caused by the tunneling: 
M{H,, t) = Mi„ + (Meq(if,) - Mi„)^r,qrt(if,)t (4.5) 

Here Min is the initial magnetization at time t = Q {i.e. after a rapid field change), 
and MQq{Hz) is the equilibrium magnetization at H^. The rate function rsqrt(-ffz) 
is proportional to the normalized distribution P{Hz) of molecules which are in 
resonance at Hz '. 

Ts^rtiH.) = (4.6) 

where is the fine width coming from the nuclear spins, Ed is the Gaussian half- 
width of P{Hz), and c is a constant of the order of unity which depends on the 
sample shape. If these simple relations are exact, then measurements of the short 
time relaxation as a function of the apphed field gives directly the distribution 
P{Hz), and allows one to measure the turmel splitting A±s which is described in 
the next section. 



4.2.2 Hole digging method to study dipolar distributions and hyperfine 
couplings 

Motivated by the Prokof'ev-Stamp theory [143J, we developed a new technique, 
which we call the hole digging method that can be used to observe the time evolution 
of molecular states in crystals of molecular clusters. It allowed us to measure the 
statistical distribution of magnetic bias fields in the Feg system that arise from the 
weak dipole fields of the clusters themselves. A hole can be 'dug' into the distribution 
by depleting the available spins at a given applied field. Our method is based on the 
simple idea that after a rapid field change, the resulting short time relaxation of the 
magnetization is directly related to the number of molecules which are in resonance 
at the given applied field. Prokof'ev and Stamp have suggested that the short time 
relaxation should follow a Vi— relaxation law (Eq. 4.5). However, the hole digging 
method should work with any short time relaxation law, for example a power law: 

M{H,,t) = Min + (Meq(F,) - Min)(r,hort(i?.)0" (4.7) 

where Fshort is a characteristic short time relaxation rate which is directly related 
to the number of molecules which are in resonance at the applied field H^, and 
< a < 1 in most cases, a = 0.5 in the Prokof'ev-Stamp theory (Eq. 4.5) and Fsqrt 
is directly proportional to P{Hz) (Eq. 4.6). The hole digging method can be divided 
into three steps (Fig. 4.12): 

1. Preparing the initial state: A well defined initial magnetization state of the 
crystal of molecular clusters can be achieved by rapidly cooling the sample 
from high down to low temperatures in a constant applied field H^. For 
zero apphed field {H^ = 0) or rather large applied fields {H^ > 1 T), one 
yields the demagnetized or saturated magnetization state of the entire crystal. 
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respectively. One can also quench the sample in a small field of few milUteslas 
yielding any possible initial magnetization Mjn. When the quench is fast (< 
1 s), the sample's magnetization does not have time to relax, either by thermal 
or quantum transitions. This procedure yields a frozen thermal equihbrium 
distribution. Whereas for slow cooling rates, the molecule spin states in the 
crystal might tend to certain dipolar ordered ground state. 

2. Modifying the initial state - hole digging: After preparing the initial state, a 
field i?dig is applied during a time t^ig^ called 'digging field and digging time', 
respectively. During the digging time and depending on -ffdig> a fraction of the 
molecular spins turmel (back and/or forth), i.e. they reverse the direction of 
magnetization. ^ 

3 . Probing the final state: Finally, a field HY°^^ is applied (Fig. 4. 1 2) to measure 
the short time relaxation from which one yields Fghort (Eq. 4.7) which is related 
to the number of spins which are still free for tunnehng after step (2). 

The entire procedure is then repeated many times but at other fields /fP™*"* yielding 

Fshort(-ffz, -ffdig, ^dig) which is related to the distribution of spins P{Hz- Hdig, tdig) 
which are still free for tunnehng after the hole digging. For t^ig = 0, this method 
maps out the initial distribution. 

4.2.3 Intermolecular dipole interaction in Fes 

We applied the hole digging method to several sample of molecular clusters and 
quantum spin glasses. The most detailed study has been done on the Fes system. 
We found the predicted ^/t relaxation (Eq. (4.5)) in experiments on fully saturated 
Fes crystals [118, 161] and on non-saturated samples [129]. Fig. 4.13 displays a 
detailed study of the dipolar distributions revealing a remarkable structure which is 
due to next nearest neighbor effects [129]. ^ These results are in good agreement 
with simulations [162, 163]. 

For a saturated initial state, the Prokof'ev-Stamp theory allows one to estimate 
the tunnel spHtting A±s- Using Eqs. (3), (9) and (12) of [143], and integration, 

we find / Fsqrt^C = ^-g^^^, where c is a constant of the order of unity which 
depends on the sample shape. With Ed = 15 mT, = 0.8 mT, c = 1 and Fsq^t from 
Fig. 4.13 in [129], we find A±io 1.2 x lO""^ K which is close to the result of 
A±io = 1.0 X 10"'^ K obtained by using a Landau-Zener method (Sect. 4. 1 .2) [47] . 



^The field sweeping rate to apply -ffdig should be fast enough to minimize the change of the initial state 
during the field sweep. 

'The peak at 0.04 T as well as the shoulder at 0.02 T and 0.04 T are originated by the clusters which have 
one nearest neighbor cluster with reversed magnetization: The peak at 0.04 T corresponds to the reversal 
of the neighboring cluster along the a crystallographic axis, which almost coincide with the easy axis of 
magnetization, while the shoulder at 0.02 T and 0.04 T are due to the clusters along b and c. 
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Fig. 4.12 Schema of the hole digging method presenting the time dependence of temperature, 
applied field and magnetization of the sample. 
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Fig. 4.13 Field dependence of the short time square-root relaxation rates Tgq^i {H,) for 
three different values of the initial magnetization -Min- According to Eq. (4.6), the curves are 
proportional to the distribution P{Hz) of magnetic energy bias due to local dipole field dis- 
tributions in the sample. Note the logarithmic scale for Fsqrt. The peaked distribution labeled 
Mia = — 0.998Ms was obtained by saturating the sample, whereas the other distributions 
were obtained by thermal anneaUng. Mi„ = — 0.870 Mg is distorted by nearest neighbor 
lattice effects. 



Whereas the hole digging method probes the longitudinal dipolar distribution (Hz 
direction), the Landau-Zener method can be used to probes the transverse dipo- 
lar distribution by measuring the tunnel splittings A around a topological quench. 
Fig. 4.14 displays such a study for the quantum transition between m = ±10, and 
m = — 10 and 9. Particular effort were made to align well the transverse field in 
direction of the hard axis. The initial magnetizations < Mjn < Mg were prepared 
by rapidly quenching the sample from 2 K in the present of a longitudinal applied 
field Hz. The quench takes approximately one second and thus the sample does not 
have time to relax, either by thermal activation or by quantum transitions, so that the 
high temperature 'thermal equilibrium' spin distribution is effectively frozen in. For 
Hz > 1 T, one gets an almost saturated magnetization state. 

The measurements of A(A'/in) show a strong dependence of the minimal tunnel 
splittings on the initial magnetization (Fig. 4.14) which demonstrates the transverse 
dipolar interaction between Feg molecular clusters being largest of Mjn = similar 
to the longitudinal dipolar interaction (Fig. 4.13). 

4.2.4 Hyperfine interaction in Fes and Mni2 

The strong influence of nuclear spins on resonant quantum tunneling in the molecular 
cluster Feg was demonstrated for the first time [144] by comparing the relaxation rate 
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Fig. 4.14 Detailed measurement of the tumiel splitting A around a topological quench for the 
quantum transition between m = — 10 and (10 — n) at = 0°. Note the strong dependence 
on the initial magnetization M-^ which demonstrates the transverse dipolar interaction between 
Feg molecular clusters [129]. 



of the standard Feg sample with two isotopic modified samples: (i) ^ Fe is replaced 
by ^'^Fe, and (ii) a fraction of is replaced by ^H. By using the hole digging method, 
we measured an intrinsic broadening which is driven by the hyperfine fields. Our 
measurements are in good agreement with numerical hyperfine calculations [144, 
136]. For T > 1.5 K, the influence of nuclear spins on the relaxation rate is less 
important, suggesting that spin-phonon coupling dominates the relaxation rate. 

Concerning Mni2 we did not find that the relaxation follows the relaxation 
law at low temperatures |164J. It is well known that the situation in this sample is 
more complicated due to the fact that there are several coexisting species of Mni2 
in any crystal, each with different relaxation times. In Ref. [164] we were able to 
isolate one faster relaxing species. The relaxation could be approximately fit to the 
^/t— relaxation law, but in fact is better fit to a power law with 0.3 < a < 0.5 
(depending on the applied field). We applied the hole digging method to this species, 
and found evidence for intrinsic line broadening below 0.3K which we suggest comes 
from nuclear spins in analogy with Peg. We also measured the relaxation of Mni2 
at higher temperature (0.04 - 5 K) and smaU fields (< 0. 1 T), and found no evidence 
for a short time ^/t— relaxation. 

4.2.5 Temperature dependence of the Landau-Zener tunneling probability 

In this section, we present studies of the temperature dependence of the Landau- 
Zener tunneling probability P yielding a deeper insight into the spin dynamics of the 
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Fcg cluster. By comparing the three isotopic samples (Sect. 4.2.4) we demonstrate 
the influence of nuclear spins on the tunneling mechanism and in particular on the 
lifetime of the first excited states. Our measurements show the need of a generalized 
Landau-Zener transition rate theory taking into account environmental effects such 
as hyperfine and spin-phonon coupling 

All measurement so far were done in the pure quantum regime (T < 0.36 K) where 
transition via excited spin levels can be neglected. We discuss now the temperature 
region of small thermal activation (T < 1 K) where we should consider transition via 
excited spin levels as well [140, 138]. 

In order to measure the temperature dependence of the tunneling probability, we 
used the Landau-Zener method as described in Sect. 4.1.2 with a phenomenological 
modification of the tunneling probabiUty P (for a negative saturated magnetization): 

P = n-wP-10,10 + Pth (4.8) 

where P_ 10, 10 is given by Eq. 4.2, n_ 10 is the Boltzmann population of them = —10 
spin level, and Pth is the overall turmeUng probability via excited spin levels. n_ lo ~ 
1 for the considered temperature T < 1 K and a negative saturated magnetization of 
the sample. 

Fig. 4. 15 displays the measured tunneling probability P for **Fe8 as a function of a 
transverse field Hx and for several temperatures. The oscillation of P are seen for 
all temperatures but the periods of oscillations decreases for increasing temperature 
(Fig. 4.17). This behavior can be explained by the giant spin model (Eq. 4.1) with 
forth order transverse terms (Sect. 4.1.3.2). Indeed, the tunnel splittings of excited 
spin levels oscillate as a function of Hx with decreasing periods (Fig. 4.16). 
Fig. 4.18 presents the tunneling probability via excited spin levels Pt^ — P — 
n_ioP-io,io- Surprisingly, the periods of Pth are temperature independent in the 
region T < 0.7 K. This suggests that only transitions via excited levels m = ±9 are 
important in this temperature regime. This statement is confirmed by the following 
estimation [166]. 

Using Eq. 4.2, typical field sweeping rates of 0.1 T/s, and tunnel splittings from 
Fig. 4.16, one easily finds that the Landau-Zener tunneling probabihty of excited 
levels are P-m,m ~ 1 for m < 10 and H 0. This means that the relaxation rates 
via excited levels are mainly governed by the lifetime of the excited levels and the 
time Tres,m during which these levels are in resonance. The later can be estimated 



" Spin-phonon interactions mainly originate from the perturbation of the crystal field by lattice vibration, 
which produce both a fluctuating local strain and a fluctuating local rotation [165, 140, 141]. It is sufficient 
to retain the lowest order terms which are quadratic with respect to spin operators. The resulting spin- 
phonon Hamiltonian contains (i) terms which commute with Sz and do not contribute to the relaxation; 
(ii) terms proportional to Sz 5+ and Sz S— ; (iii) terms proportional to and . Thus, the spin-phonon 
interaction has matrix elements between states with quantum numbers m and m' if |m — m'\ = 1 or 
2 [140]. 
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Fig. 4.15 Transverse field dependence of the tunneling probability P at several temperatures, 
and the ground state tunneling probability P_ 10,10 measured at T = 0.05 K and for "'Fes. 
The field sweeping rate was 0.14 T/s. The dotted Unes indicate the minima of P_io^io. 




Fig. 4.16 Calculated tunnel splitting A^^jn' as a function of the transverse field for 
quantum transition between m = ±10, ±9 and ±8 (Sect. 4.1.3.2). 
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Fig. 4.17 Temperature dependence of the field between minima of the tunneling probability 
P in Fig. 4.15. a, b, and c are defined in the inset. The dotted line labeled with a', b', and c' 
where take from Pth of Fig. 4.18; see also [126]. 



10-5 

-0.4 " # ' m " .0^:; ' -iM 

Fig. 4.18 Transverse field dependence of Pih which is the difference between the measured 
tunnel probability P and the ground state tunnel probability Tt-io^'-io.lO measured at T = 
0.05K (see Fig. 4. 15. The field sweeping rate was 0.14 T/s. The long dotted lines indicate the 
minima of Pth whereas the short dotted Unes indicate the minima of P-io,io- 
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by 

guBmnodHz/dt 

The probability for a spin to pass into the excited level m can be estimated by 
T^^e~^^°-'"^''^'^, where -Eio.m is the energy gap between the levels 10 and m, and 
Tm is the lifetime of the excited level m. One gets: 



Pth ~ Tres,m ^-Eio,^/ksT ^ ^-m,m ^- 



Eio,m/feBT 



(4.10) 

Note that this estimation neglects higher excited levels with |m| < 8. ^' Fig. 4.19 
displays the measured Pth for the three isotopic Feg samples. For 0.4 K < T < 1 K 
we fitted Eq. 4.10 to the data leaving only the level lifetimes ro and rs as adjustable 
parameters. All other parameters are calculated using the parameters in Sect. 4.1.3.2. 
We obtain rg = 1.0,0.5, and 0.3 x lO'^^s, and ts = 0.7,0.5, and 0.4 x lO-'^s 
for ^Fes, '*Fe8, and ^^Fes, respectively. These results evidence that only the first 
excited level have to be considered for 0.4 K < T < 0.7 K. Indeed, the second 
term of the summation in Eq. 4.10 is negligible in this temperature interval. It is 
interesting to note that this finding is in contrast to hysteresis loop measurements on 
Mni2 [120, 167] which suggested an abrupt transition between thermal assisted and 
pure quantum tunneling [168]. Furthermore, our result shows clearly the influence 
of nuclear spins which seem to decrease the level lifetimes Tm, i-e. to increase 
dissipative effects. 

The nuclear magnetic moment and not the mass of the nuclei seems to have the major 
effect on the dynamics of the magnetization. In fact the mass is increased in both 
isotopically modified samples whereas the effect on the relaxation rate is opposite. 
On the other hand ac susceptibihty measurements at T > 1 .5 K showed no clear 
difference between the three samples suggesting that above this temperature, where 
the relaxation is predominately due to spin-phonon coupUng [140, 141], the role 
of the nuclear spins is less important. Although the increased mass of the isotopes 
changes the spin-phonon coupling, this effect seems to be small. 
We can also exclude that the change of mass for the three isotopic samples has 
induced a significant change in the magnetic anisotropy of the clusters. In fact the 
measurements below T < 0.35 K, where spin-phonon coupling is negligible, have 
shown that (i) relative positions of the resonances as a function of the longitudinal field 
Hz are unchanged and (ii) all three samples have the same period of oscillation 



^'The probability of phonon induced transitions with | Am| > 2 are very small [139, 140, 138], Also the 
Boltzmann factor e-^^o.m/knT ygry small for m < 8 and T < 1 K. 

'^We observed a small shift of the resonances of the order of magnitude of 1 mT, positive for ^'^Fes and 
negative for '^Feg (Minit = —Ms). This can also be attributed to the modified hyperfine fields. However, 
a quantitative measurement is complicated by the fact that it is impossible to have two crystal with exactly 
the same shape, i.e. the same internal fields. 
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Fig. 4.19 Temperature dependencies of Pth for Hx = for three Peg samples. The field 
sweeping rate was 0.14 T/s. The dotted lines are fits of the data using Eq. 4.10. 



of A as a function of the transverse field [47], a period which is very sensitive to 
any change of the anisotropy constants. 

4.2.6 Conclusion on molecular magnets 

In conclusion, we presented detailed measurements which demonstrated that molec- 
ular magnets offer an unique opportunity to explore the quantum dynamics of a large 
but finite spin. We focused our discussion on the Fes molecular magnet because 
it is the first system where studies in the pure quantum regime were possible. The 
tunneling in this system is remarkable because it does not show up at lowest order of 
perturbation theory. 

What remains still debated is the possibility of observing quantum coherence between 
states of opposite magnetization. Dipole-dipole and hyperfine interactions are source 
of decoherence. In other words when a spin has tunneled through the barrier, it 
experiences a huge modification of its environment (hyperfine and dipolar) which 
prohibits the back tunneling. Prokof ' ev and Stamp suggested three possible strategies 
to suppress the decoherence [172]. (i) Choosing a system were the NMR frequencies 
far exceed the tunnel frequencies making any coupling impossible, (ii) Isotopically 

'^Molecules with small spin have also been studied. For example, time resolved magnetization measure- 
ments were performed on a spin 1/2 molecular complex, so-called Vis [169]. Despite the absence of 
a barrier, magnetic hysteresis is observed over a time scale of several seconds. A detailed analysis in 
terms of a dissipative two-level model has been given, in which fluctuations and splittings are of the same 
energy. Spin-phonon coupling leads to long relaxation times and to a particular 'butterfly' hysteresis 
loop [170, 171]. 
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purify the sample to remove all nuclear spins, (iii) Applying a transverse field to 
increase the tunnel rate to frequencies much larger than hyperfine field fluctuations. 
All three strategies are difficult to realize. However, some authors tried to realize 
the last one by performing EPR experiments in the presence of a large transverse 
field [173]. Absorption of radio-frequency electromagnetic fields were observed 
which might be due to induced transitions near the tunnel splitting. However, no 
experiments showed the oscillatory behavior in the time domain which might be 
evidenced by a spin echo type of experiment. 

Concerning the perspectives of the field of single molecule magnets we expect that 
chemistry is going to play a major role through the synthesis of novel larger spin 
clusters with strong anisotropy. We want to stress that there are already many other 
molecular magnets (see for instance [121, 122, 123, 124]) which are possible model 
systems. We beUeve that more sophisticated theories are needed which describe 
the dephasing effects of the environment onto the quantum system. These studies 
are important for sing the quantum character of molecular clusters for applications 
like 'Quantum Computers' . The first implementation of Grover's algorithm with 
molecular magnets has been proposed [174] 

4.3 QUANTUM TUNNELEVG OF MAGNETIZATION IN ESfDIVIDUAL 
SINGLE-DOMAIN NANOPARTICLES 

The following sections focuses on Magnetic Quantum Tunneling (MQT) studied in 
individual nanoparticles or nanowires where the complications due to distributions 
of particle size, shape, and so on are avoided. The experimental evidence for MQT in 
a single-domain particle or in assemblies of particles is still a controversial subject. 
We shall therefore concentrate on the necessary experimental conditions for MQT 
and review some experimental results which suggest that quantum effects might even 
be important in nanoparticles with S — 10^ or larger. We start by reviewing some 
important predictions concerning MQT in a single-domain particle. 

4.3.1 Magnetic Quantum TUnneling in nanoparticles 

On the theoretical side, it has been shown that in small magnetic particles, a large 
number of spins coupled by strong exchange interaction, can tunnel through the 
energy barrier created by magnetic anisotropy. It has been proposed that there is a 
characteristic crossover temperature Tc below which the escape of the magnetization 
from a metastable state is dominated by quantum barrier transitions, rather than by 
thermal over barrier activation. Above Tc the escape rate is given by thermal over 
barrier activation (Sect. 3). 

In order to compare experiment with theory, predictions of the crossover temperature 
Tc and the escape rate Fqt in the quantum regime are relevant. Both variables should 
be expressed as a function of parameters which can be changed experimentally. 
Typical parameters are the number of spins S, effective anisotropy constants, appUed 
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Fig. 4.20 Normalized crossover temperature as given by (4. 11), and for several values of 



field strength and direction, coupling to the environment (dissipation) etc. Many 
theoretical papers have been published during the last few years [108]. We discuss 
here a result specially adapted for single particle measurements, which concerns the 
field dependence of the crossover temperature T^. 

The crossover temperature Tc can be defined as the temperature where the quantum 
switching rate equals the thermal one. The case of a magnetic particle, as a function 
of the applied field direction, has been considered by several authors [175, 176, 177]. 
We have chosen the result for a particle with biaxial anisotropy as the effective 
anisotropy of most particles can be approximately described by strong uniaxial and 
weak transverse anisotropy. The result due to Kim can be written in the following 
form [177]: 



where ^,oH\\ = K^^/M^ and iioH±_ = K±/Ms are the parallel and transverse 
anisotropy fields given in Tesla, K^^ and K± are the parallel and transverse anisotropy 
constants of the biaxial anisotropy, 9 is the angle between the easy axis of magneti- 
zation and the direction of the applied field, ande = {1 — H/H^^). Equation (4.11) 
is valid for any ratio a = H^/H\y The proportionality coefficient of (4.11) is of 
the order of unity (Tc is in units of Kelvin) and depends on the approach used for 
calculation [177]. Equation (4.11) is plotted in Fig. 4.20 for several values of the 
ratio n. It is valid in the range \/e < 6 < tt /2 — ^fe. 

The most interesting feature which may be drawn from (4.11) is that the crossover 
temperature is tunable using the external field strength and direction (Fig. 4.20) 
because the tunneUng probability is increased by the transverse component of the ap- 



the ratio a = H±/H\\ 




(4.11) 
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plied field. Although at high transverse fields, decreases again due to a broadening 
of the anisotropy barrier. Therefore, quantum tunneling experiments should always 
include studies of angular dependencies. When the effective magnetic anisotropy of 
the particle is known, MQT theories give clear predictions with no fitting parameters. 
MQT could also be studied as a function of the effective magnetic anisotropy. In 
practice, it is well known for single particle measurements that each particle is some- 
what different. Therefore, the effective magnetic anisotropy has to be determined for 
each particle (Sect. 2.1.1). 

Finally, it is important to note that most of the MQT theories neglect damping 
mechanisms. We discussed the case of ohmic damping in Sect. 3.1 [90] which is the 
simplest form of damping. More complicated damping mechanisms might play an 
important role. We expect more theoretical work on this in future. 

4.3.2 Magnetization measurements of individual single-domain nanoparticles 
and wires at very low temperatures 

In order to avoid the complications due to distributions of particle size, shape, and so 
on, some groups have tried to study the temperature and field dependence of magne- 
tization reversal of individual magnetic particles or wires. Most of the recent studies 
were done using Magnetic Force Microscopy at room temperature. Low temperature 
investigations were mainly performed via resistance measurements (Sect. 1.1). 
The first magnetization measurements of individual single-domain nanoparticles at 
low temperature (0.1 K < T < 6 K) were presented by Wernsdorfer et al. [22]. 
The detector (a Nb micro-bridge-DC-SQUID, see Sect. 1 .2) and the particles studied 
(ellipses with axes between 50 and 1000 nm and thickness between 5 and 50 nm) were 
fabricated using electron beam hthography. Electrodeposited wires (with diameters 
ranging from 40 to 100 nm and lengths up to 5000 nm) were also studied [35, 
74]. Waiting time and switching field measurements (Sect. 3.2) showed that the 
magnetization reversal of these particles and wires results from a single thermally 
activated domain wall nucleation, followed by a fast wall propagation reversing 
the particle's magnetization. For nanocrystalline Co particles of about 50 nm and 
below 1 K, a flattening of the temperature dependence of the mean switching field 
was observed which could not be explained by thermal activation. These results 
were discussed in the context of MQT. However, the width of the switching field 
distribution and the probabihty of switching are in disagreement with such a model 
because nucleation is very sensitive to factors Uke surface defects, surface oxidation 
and perhaps nuclear spins. The fine structure of pre-reversal magnetization states is 
then governed by a multivalley energy landscape (in a few cases distinct magnetization 
reversal paths were effectively observed [59]) and the dynamics of reversal occurs 
via a complex path in configuration space. 

Coppinger et al. [23] used telegraph noise spectroscopy to investigate two-level 
fluctuations (TLF) observed in the conductance of a sample containing self-assembled 
ErAs quantum wires and dots in a semi-insulating GaAs matrix. They showed that 
the TLF could be related to two possible magnetic states of a ErAs cluster and that 
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the energy difference between the two states was a linear function of the magnetic 
field. They deduced that the ErAs cluster should contain a few tens of Er atoms. 
At temperatures between 0.35 K and 1 K, the associated switching rates of the TLF 
were thermally activated, whilst below 0.35 K the switching rate became temperature 
independent. Turmeling of the magnetization was proposed in order to explain the 
observed behavior. 

Some open questions remain: what is the object which is really probed by TLF? If this 
is a single ErAs particle, as assumed by the authors, the switching probability should 
be an exponential function of time. The pre-exponential factor Tq~^ (sometimes called 
attempt frequency) was found to lie between 10^ and 10^ s~ ^ whereas expected values 
are between 10® and 10^^ s~^. Why must one apply fields of about 2 T in order to 
measure two-level fluctuations which should be expected near zero field? What is 
the influence of the measurement technique on the sample? 

By measuring the electrical resistance of isolated Ni wires with diameters between 20 
and 40 nm, Hong and Giordano studied the motion of magnetic domain walls [24]. 
Because of surface roughness and oxidation, the domain walls of a single wire are 
trapped at pinning centers. The pinning barrier decreases with an increase in the 
magnetic field. When the barrier is sufficiently smaU, thermaUy activated escape of 
the wall occurs. This is a stochastic process which can be characterized by a switching 
(depinning) field distribution. A flattening of the temperature dependence of the mean 
switching field and a saturation of the width of the switching field distribution (rms. 
deviation a) were observed below about 5 K. The authors proposed that a domain 
wall escapes from its pinning site by thermal activation at high temperatures and by 
quantum tunneling below Tc 5 K. 

These measurements pose several questions: what is the origin of the pinning center 
which may be related to surface roughness, impurities, oxidation and so on? The 
sweeping rate dependence of the depinning field, as well as the depinning probability, 
could not be measured even in the thermally activated regime. Therefore, it was not 
possible to check the validity of the Neel-Brown model [9, 10, 86, 87, 88] or to 
compare measured and predicted rms. deviations a. Finally, a crossover temperature 
Tc of about 5 K is three orders of magnitude higher than Tc predicted by current 
theories. 

Later, Wemsdorfer et al. published results obtained on nanoparticles synthesized 
by arc discharge, with dimensions between 10 and 30 nm [36]. These particles 
were single crystalUne, and the surface roughness was about two atomic layers. 
Their measurements showed for the first time that the magnetization reversal of a 
ferromagnetic nanoparticle of good quality can be described by thermal activation 
over a single-energy barrier as proposed by Neel and Brown [9, 10, 86, 87, 88] (see 
Sect. 3.3). The activation volume, which is the volume of magnetization overcoming 
the barrier, was very close to the particle volume, predicted for magnetization reversal 
by uniform rotation. No quantum effects were found down to 0.2 K. This was not 
surprising because the predicted cross-over temperature is Tc 0.02 K. The results 
of Wernsdorfer et al. constitute the preconditions for the experimental observation 
of MQT of magnetization on a single particle. 
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Fig. 4.21 Angular dependence of the crossover temperature Tc for a 10 nm 
BaFe10.4Coo.8Tio.8O19 particle with S « 10^. The lines are given by (4.11) for differ- 
ent values of the ratio a = H±/H\\. The experimental data are normalized by rc(45°) = 
0.31 K 



Just as the results obtained with Co nanoparticles [36], a quantitative agreement with 
the Neel-Brown model of magnetization reversal was found on BaFei2-2xCOa;Ti2:Oi9 
nanoparticles(0 < a; < 1) [37], which we will call BaFeO, in the size range of 10 - 
20 nm. However, strong deviations from this model were evidenced for the smallest 
particles containing about lO'^/ie and for temperatures below 0.4 K. These devia- 
tions are in good agreement with the theory of macroscopic quantum tunneUng of 
magnetization. Indeed, the measured angular dependence of Tc{6) is in excellent 
agreement with the prediction given by (4.1 1) (Fig. 4.21). The normalization value 
T(,(45°) = 0.31 K compares well with the theoretical value of about 0.2 K. 
Although the above measurements are in good agreement with MQT theory, we 
should not forget that MQT is based on several strong assumptions. Among them, 
there is the assumption of a giant spin, i.e., all magnetic moments in the particle are 
rigidly coupled together by strong exchange interaction. This approximation might 
be good in the temperature range where thermal activation is dominant but is it not 
yet clear if this can be made for very low energy barriers. Future measurements 
might give us the answer. 

4.3.3 Quantization of the magnetization 

In order to give a definite proof that MQT can occur in a magnetic nanoparticle we 
propose to surge for the energy level quantization of its collective spin state. This was 
recently evidenced in molecular cluster Uke Fe8 having a collective spin state S' = 10 
(Sect. 4.1). In the case of the BaFeO particles with S « 10^ [37], the field separation 
associated with level quantization is rather small: AH w Hg,/2S 0.002 mT where 
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ifa is the anisotropy field. However, for a 3 nm Co cluster with S ~ 10"^ the field 
separation AH = H^/2S ~ 0.2 mT might be large enough to be measurable. 
Fig. 4.22 displays schematically the field values of resonances between quantum 
states of S. When the appUed field is ramped in a certain direction, the resonance 

might occur for fields i?res ~ ^ ' with n = 1 , 2, 3 , ^ is the angle between 

the applied field and the easy axis of magnetization. For large spins S, tunneling 
might be observable only for fields which are close to the classical switching field 
(Fig. 4.22). The resonance fields could be evidenced by measuring switching field 
distributions (inset of Fig. 4.22) as a function of the angle 9. 

Such a study is presented in Fig. 4.23 for a 3 nm Fe cluster with S « 800. The 
estimated field separation Aif = Hg_/2S is about 0.1 mT whereas the width of 
the switching field distribution is about ten times larger. We observed sometimes a 
small periodic fine structure which is close to the expected AH. However this fine 
structure always disappeared when averaging over more measurements. A possible 
origin might be hyperfine couplings which broaden the energy levels (Sect. 4.2.4) 
leading to a complete overlap of adjacent energy levels. It is also important to mention 
that the switching field distributions were temperature independent for 0.04 K < T < 
0.2 K. 

In some cases, we observed huge variations of the switching field (Figs. 3.6 and 3.7) 
which might be due to exchange bias of frustrated spin configurations. However, 
quantum effects are not completely excluded. 

Future measurements should focus on the level quantization of collective spin states 
ofS' = 102. 
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Fig. 4.22 Schematic view of the resonance fields of a giant spin S. The continuous line is the 
classical switching fields of Stoner-WoUfarth (Sect. 2.1). The inset presents schematically a 
switching field histogram with AH' ss ^ 




Fig. 4.23 Angular dependence of the switching field of a 3 nm Fe cluster with S Ri 800. The 
inset presents detailed measurements near the easy axis of magnetization. 



Summary and conclusion 

Nanometer-sized magnetic particles have generated continuous interest as the study of 
their properties has proved to be scientifically and technologically very challenging. 
In the last few years, new fabrication techniques have led to the possibility of making 
small objects with the required structural and chemical qualities. In order to study 
these objects new local measuring techniques were developed such as magnetic force 
microscopy, magnetometry based on micro-Hall probes or micro-SQUIDs. This led 
to a new understanding of the magnetic behavior of nanoparticles. 
In this article we reviewed the most important theories and experimental results con- 
cerning the magnetization reversal of single-domain particles and clusters. Special 
emphasis is laid on single particle measurements avoiding complications due to dis- 
tributions of particle size, shape etc. Measurements on particle assembUes has been 
reviewed in [8] . We mainly discuss the low temperature regime in order to avoid spin 
excitations. 

Sect. 1 reviews briefly the commonly used measuring techniques. Among them, 
electrical transport measurements. Hall probes and micro-SQUID techniques seem 
to be the most convenient techniques for low temperature measurements. 
Sect. 2 discusses the mechanisms of magnetization reversal in single domain particles 
at zero Kelvin. For extremely small particles, the magnetization should reverse by 
uniform rotation of magnetization (Sec. 2.1). For somewhat larger particles, non 
uniform reversal mode occur like the curling mode. (Sec. 2.2). For even larger 
particles, magnetization reversal occurs via a domain wall nucleation process starting 
in a rather small volume of the particle (Sec. 2.2.3). 

The influence of temperature on the magnetization reversal is reported in Sect. 3. We 
discuss in detail the Neel, Brown and Coffey's theory of magnetization reversal by 
thermal activation (Sect. 3). 

Finally, Sect. 4 shows that for very small systems or very low temperature, mag- 
netization can reverse via quantum tunneling. The boundary between classical and 
quantum physics has become a very attractive field of research. This section dis- 
cusses detailed measurements which demonstrated that molecular magnets offer an 
unique opportunity to explore the quantum dynamics of a large but finite spin. The 
discussion is focused on the Feg molecular magnet with 5 = 10 because it is the first 
system where studies in the pure quantum regime were possible. We showed that the 
understanding of the environmental decoherence is one of the most important issues, 
in particular for future applications of quantum devices. We then discussed tunneling 
in nanoparticles and showed how one might give a definite proof of their quantum 
character at low temperature. 
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In conclusion, the understanding of the magnetization reversal in nano-structures 
requires the knowledge of many physical phenomena, and are therefore particularly 
interesting for the development of new fundamental theories of magnetism and in 
modeling new magnetic materials for permanent magnets or high density record- 
ing. Using the quantum character of nano-structures for applications like 'Quantum 
Computers' wiU be one of the major concerns of the next decades. 
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